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Abstract 

A search for flavour-changing neutral current decays of a top quark to an up-type quark {q = 
u, c) and the Standard Model Higgs boson, where the Higgs boson decays to bb, is presented. 
The analysis searches for top quark pair events in which one top quark decays to Wb, with the 
W boson decaying leptonically, and the other top quark decays to Hq. The search is based 
on pp collisions at = 8 TeV recorded in 2012 with the ATLAS detector at the CERN 
Large Hadron Collider and uses an integrated luminosity of 20.3 fb”^ Data are analysed in 
the lepton-plus-jets final state, characterised by an isolated electron or muon and at least four 
jets. The search exploits the high multiplicity of ^-quark jets characteristic of signal events, 
and employs a likelihood discriminant that uses the kinematic differences between the signal 
and the background, which is dominated by tt —> WbWb decays. No significant excess of 
events above the background expectation is found, and observed (expected) 95% CL upper 
limits of 0.56% (0.42%) and 0.61% (0.64%) are derived for the t He and t —> Hu 
branching ratios respectively. The combination of this search with other ATLAS searches in 
the H ^ yy and H WW*,tt decay modes significantly improves the sensitivity, yielding 
observed (expected) 95% CL upper limits on the t He and t —> Hu branching ratios 
of 0.46% (0.25%) and 0.45% (0.29%) respectively. The corresponding combined observed 
(expected) upper limits on the \AtcH\ and \AtuH\ couplings are 0.13 (0.10) and 0.13 (0.10) 
respectively. These are the most restrictive direct bounds on tqH interactions measured so 
far. 
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1. Introduction 

Following the observation of a Higgs boson by the ATLAS and CMS collaborations [1, 2], a comprehens¬ 
ive programme of measurements of its properties is underway looking for deviations from the Standard 
Model (SM) predictions. An interesting possibility is the presence of flavour-changing neutral current 
(FCNC) interactions between the Higgs boson, the top quark, and a u- or c-quark, tqH {q = u, c). Since 
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the Higgs boson is lighter than the top quark, with a measured mass mn - 125.09 + 0.24 GeV [3], such 
interactions would manifest themselves as FCNC top quark decays, t Hq. In the SM, such decays are 
extremely suppressed relative to the dominant t —> Wb decay mode, since tqH interactions are forbidden 
at the tree level and even suppressed at higher-orders in the perturbative expansion due to the Glashow- 
Iliopoulos-Maiani (GIM) mechanism [4]. As a result, the SM predictions for the t Hq branching 
ratios are exceedingly small: BR(t ^ Hu) ~ 10“^^ and BR(t —> He) ~ 10“^^ [5-8]. On the other hand, 
large enhancements in these branching ratios are possible in some beyond-SM scenarios, where the GIM 
suppression can be relaxed and/or new particles can contribute to the loops, yielding effective couplings 
orders of magnitude larger than those of the SM. Examples include quark-singlet models [9], two-Higgs- 
doublet models (2HDM) of type I, with explicit flavour conservation, and of type II, such as the minimal 
supersymmetric SM (MSSM) [10-12], or supersymmetric models with R-parity violation [13]. In those 
scenarios, typical branching ratios can be as high as BR(f ^ Hq) ~ 10“^. An even larger branching ratio 
of BR(f ^ He) ~ 10~^ can be reached in 2HDM without explicit flavour conservation (type III), since a 
tree-level FCNC coupling is not forbidden by any symmetry [14-16]. While other FCNC top couplings, 
tqy, tqZ, tqg, are also enhanced relative to the SM prediction in those scenarios beyond the SM, the largest 
enhancements are typically for the tqH couplings, and in particular the tcH coupling. See ref. [7] for a 
review. 

Searches for t Hq decays have been performed by the ATFAS and CMS collaborations, taking ad¬ 
vantage of the large samples of tt events collected during Run 1 of the FHC. In these searches, one 
of the top quarks is required to decay into Wb, while the other top quark decays into Hq, yielding 
tt — > WbHq} Assuming SM decays for the Higgs boson and tue = 125 GeV, the most sensitive single¬ 
channel searches have been performed in the H ^ yy decay mode which, despite the tiny branching 
ratio of BR(// ^ yy) - 0.2%, is characterised by very small background and excellent diphoton mass 
resolution. The resulting observed (expected) 95% confidence level (CF) upper limits on BR(f ^ Hq) 
are 0.79% (0.51%) and 0.69% (0.81%), respectively from the ATFAS [17] and CMS [18] collaborations. 
These searches are insensitive to the difference between t Hu and t He, and thus the above limits can 
be interpreted as applying to the sum BR(t ^ Hu)+BR(t He). The CMS Collaboration has also rein¬ 
terpreted searches in multilepton (three or four leptons) final states [18] in the context of tt WbHq with 
H WW* ,TT, resulting in an observed (expected) upper limit of BR(f —> He) < 1.28% (1.17%) at the 
95% CF. Multilepton searches are able to exploit a significantly larger branching ratio for the Higgs boson 
decay compared to the H ^ yy decay mode, and are also characterised by relatively small backgrounds. 
However, in general they do not have good mass resolution,^ so any excess would be hard to interpret as 
originating from t —> Hq decays. The combination of CMS searches in diphoton and multilepton (three 
or four leptons) final states yields an observed (expected) upper limit of BR(f —> He) < 0.56% (0.65%) 
at the 95%CF[18]. 

Upper limits on the branching ratios BR(t ^ Hq) (q = u, e) can be translated to upper limits on the 
non-flavour-diagonal Yukawa couplings AtqH appearing in the following Fagrangian: 


-Cfcnc = ^icHtHe + AfuHtHu + h.e. (1) 

The branching ratio BR(f —> Hq) is estimated as the ratio of its partial width [8] to the SM t —> Wb partial 
width [19], which is assumed to be dominant. Both predicted partial widths include next-to-leading-order 


* In the following WbHq is used to denote both W*bHq and its charge conjugate, HqW b. Similarly, WbWb is used to denote 
W+bW-b. 

^ An exception is the H —> ZZ* —> (£, I' = e,p) decay mode, which has a very small branching ratio and thus is not 

promising for this search. 
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(NLO) QCD corrections. Using the expression derived in ref. [17], the coupling \AtqH\ can be extracted as 
\AtqH\ = (1-92 ± 0.02) VBR(t 4 Hq). 

The results presented in this paper fill a gap in the current programme of searches for t —> Hq decays at 
the LHC by considering the dominant decay mode H bb, which has BR(// ^ bb) ^ 58%. This search 
is focused on the tt WbHq {q - u, c) process, with W ^ bv {£ = e,)!, t) and H bb, resulting in a 
lepton-plus-jets final sfafe wifh high f7-jef mulfiplicify, which can be effecfively exploited fo suppress fhe 
overwhelming tt background. Early sfudies of fhe prospecfs for fhis search af fhe LHC were performed in 
ref. [20]. Only evenfs wifh an elecfron or muon, including fhose produced via lepfonically decaying fans, 
are considered. The lepfon-plus-Jefs final sfafe also allows fhe kinematic reconsfrucfion of fhe final sfafe 
and in particular fhe dijef invarianf mass specfrum from the H ^ bb decay, providing additional handles 
that would help in detecting tt —> WbHq events. Most of this paper is devoted to the discussion of this 
particular search, for which background estimation techniques, systematic uncertainties and statistical 
treatment closely follow those used in recent ATLAS searches using the same final-state signature [21, 
22] . This paper also includes a reinterpretation of the ATLAS search for ttH associated production, with 
H WW* ,ZZ* ,TT, resulting in multilepton final sfafes [23]. This reinferprelafion only considers fhe 
final sfafes with a significant expected contribution from tt WbHq, H WW*,tt signal, namely two 
same-charge leptons with and without an identified hadronic fau lepfon and fhree lepfons. A combination 
of fhe fhree ATLAS searches for tt WbHq, probing fhe H bb, H —> WW*, tt, and H ^ yy decay 
modes, is performed and bounds are sef on BR(t ^ He) and BR(t ^ Hu), as well as on fhe corresponding 
non-flavour-diagonal Yukawa couplings. 

This paper is organised as follows. A brief descripfion of fhe ATLAS defeefor is provided in seefion 2. 
Subsequenf seefions are devoted fo a defailed discussion of fhe tt WbHq, H ^ bb search, covering fhe 
objeef reconsfrucfion (seefion 3), fhe dafa sample and even! preselecfion (section 4), fhe modelling of fhe 
backgrounds and fhe signal (seefion 5), fhe analysis sfrafegy (section 6 ), and fhe sysfemafic uncerfainfies 
(seefion 7). Seefion 8 provides a discussion of fhe sfafisfical methods used. Section 9 presents the results 
obtained by the three individual ATLAS searches as well as their combination. Linally, the conclusions 
are given in section 10 . 


2. ATLAS detector 

The ATLAS detector [24] consists of the following main subsystems: an inner tracking system, elec¬ 
tromagnetic and hadronic calorimeters, and a muon spectrometer. The inner detector provides tracking 
information from silicon pixel and microstrip detectors in the pseudorapidity^ range I 77 I < 2.5 and from 
a straw-tube transition radiation tracker covering \q\ < 2.0, all immersed in a 2 T axial magnetic field 
provided by a superconduefing solenoid. The elecfromagnefic (EM) sampling calorimeter uses lead as 
fhe absorber material and liquid-argon (LAr) as fhe aefive medium, and is divided info barrel (|? 7 | < 1.475) 
and end-cap (1.375 < I 77 I < 3.2) regions. Hadron calorimefry is also based on the sampling technique, 
with either scintillator tiles or LAr as the active medium, and with steel, copper, or tungsten as the ab¬ 
sorber material. The calorimeters cover \q\ < 4.9. The muon spectrometer measures the deflection of 
muons with \q\ < 2.7 using multiple layers of high-precision tracking chambers located in a toroidal 


^ ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector 
and the z-axis coinciding with the axis of the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the 
^-axis points upward. Cylindrical coordinates (r,0) are used in the transverse plane, 0 being the azimuthal angle around the 
beam pipe. The pseudorapidity is defined in terms of the polar angle 6 as 77 = - In tan(6/2). 
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field of approximately 0.5 T and 1 T in the central and end-cap regions of ATLAS, respectively. The 
muon spectrometer is also instrumented with separate trigger chambers covering I 77 I < 2.4. A three-level 
trigger system [25] is used to select interesting events. The first-level trigger is implemented in custom 
electronics and uses a subset of detector information to reduce the event rate to at most 75 kHz. This 
is followed by two software-based trigger levels exploiting the full detector information and yielding a 
typical recorded event rate of 400 Hz during 2012. 


3. Object reconstruction 

Electron candidates [26] are reconstructed from energy clusters in the EM calorimeter that are matched 
to reconstructed tracks in the inner detector. Electron clusters are required to have a transverse energy Ex 
greater than 25 GeV and iT/dusterl < 2.47, excluding the transition region 1.37 < |77ciusterl < 1-52 between 
sections of the EM calorimeter. The longitudinal impact parameter of the electron track with respect to 
the event’s primary vertex (see section 4), zo, is required to be less than 2 mm. Electrons are required 
to satisfy “tight” quality requirements [26] based on calorimeter, tracking and combined variables that 
provide good separation between prompt electrons and jets. To reduce the background from non-prompt 
electrons resulting from semileptonic decays of b- or c-hadrons, and from jets with a high fraction of their 
energy deposited in the EM calorimeter, electron candidates must also satisfy calorimeter- and track-based 
isolation requirements. The calorimeter isolation variable is based on the energy sum of cells within a 
cone of size AT? = '^/(A^)^ - 1 - (A?/)^ = 0.2 around the direction of each electron candidate, and an rj- 
dependent requirement is made, giving an average efficiency of 90% across 77 for prompt electrons from 
Z boson decays. This energy sum excludes cells associated with the electron cluster and is corrected 
for leakage from the electron cluster itself as well as for energy deposits from additional pp interactions 
within the same bunch crossing (“pileup”). A further 90%-efficient isolation requirement is made on the 
track transverse momentum (px) sum around the electron (excluding the electron track itself) in a cone 
of size A/? = 0.3. 

Muon candidates [27, 28] are reconstructed from track segments in the various layers of the muon spec¬ 
trometer that are matched with tracks found in the inner detector. The final candidates are refitted using 
the complete track information from both detector systems and are required to have px > 25 GeV and 
Ipl < 2.5. The longitudinal impact parameter of the muon track with respect to the primary vertex, zo, is 
required to be less than 2 mm. Muons are required to satisfy a px-dependent track-based isolation require¬ 
ment: the scalar sum of the px of the tracks within a cone of variable size AR = 10 GeV/p(|. around the 
muon (excluding the muon track itself) must be less than 5% of the muon px (p^). This requirement has 
good signal efficiency and background rejection even under high-pileup conditions, as well as in boosted 
configurations where the muon is close to a jet. Eor muons from W boson decays in simulated tt events, 
the average efficiency of the isolation requirement is about 95%. 

Jets are reconstructed with the anti-kf algorithm [29-31] with a radius parameter R = 0.4, using cali¬ 
brated topological clusters [32, 33] built from energy deposits in the calorimeters. Prior to jet finding, 
a local cluster calibration scheme [34] is applied to correct the topological cluster energies for the non¬ 
compensating response of the calorimeter, as well as for the energy lost in dead material and via out-of- 
cluster leakage. The corrections are obtained from simulations of charged and neutral particles. After 
energy calibration [35], jets are required to have px > 25 GeV and \r]\ < 2.5. To reduce the contamination 
due to jets originating from pileup interactions, a requirement on the absolute value of the jet vertex 
fraction (JVE) variable above 0.5 is applied to jets with px < 50 GeV and \r]\ < 2.4. This requirement 
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ensures that at least 50% of the scalar sum of the p^; of the tracks with pi > \ GeV associated with a jet 
comes from tracks originating from the primary vertex. During jet reconstruction, no distinction is made 
between identified electrons and jet energy deposits. Therefore, if any of the jets lie within AT? = 0.2 
of a selected electron, the closest jet is discarded in order to avoid double-counting of electrons as jets. 
Finally, any electron or muon within AT? = 0.4 of a selected jet is discarded. 

Jets containing Tr-hadrons are identified (Tr-fagged) via an algorifhm [36] fhaf uses mulfivariafe fechniques 
fo combine informafion from fhe impacf paramefers of displaced fracks as well as topological properties of 
secondary and ferfiary decay verfices reconstrucfed wifhin fhe jef. For each jef, a value for fhe mulfivariafe 
Tr-fagging discriminanf is calculafed. The jef is considered Tr-fagged if fhis value is above a given fhreshold. 
The fhreshold used in fhis search corresponds fo 70% efficiency fo fag a T?-quark jef, wifh a lighf-jef^ 
rejection factor of ~130 and a charm-jef rejection factor of 5, as determined for jefs wifh pj > 20 GeV 
and \r]\ < 2.5 in simulated tt evenfs. 

The missing fransverse momenfum is consfrucfed [37] from fhe vecfor sum of all calorimefer 

energy deposifs confained in fopological clusters. All fopological clusfer energies are correcfed using fhe 
local clusfer calibrafion scheme discussed previously in fhe confexf of fhe jef energy calibrafion. Those 
fopological clusfers associafed wifh a high-px objecf (e.g. jef or elecfron) are furfher calibrated using fheir 
respecfive energy correcfions. In addifion, confribufions from fhe px of selecfed muons are included in 
fhe calculation of 


4. Data sample and event preselection 

This search is based on pp collision dafa af = 8 TeV collecfed by fhe ATLAS experimenf befween 
April and December 2012. Only evenfs recorded wifh a single-election or single-muon frigger under 
sfable beam conditions and for which all defector subsystems were operafional are considered. The cor¬ 
responding integrated luminosify is 20.3+0.6 fb“^ [38]. Single-lepfon tiiggers wifh differenf px fhresholds 
are combined in a logical OR in order fo increase fhe overall efficiency. The px fhresholds are 24 or 
60 GeV for fhe elecfron tiiggers and 24 or 36 GeV for fhe muon tiiggers. The tiiggers wifh fhe lower px 
fhreshold include isolafion requiremenfs on fhe candidate lepfon, resulfing in inefficiencies af high px thaf 
are recovered by fhe tiiggers wifh higher px fhreshold. 

Evenfs satisfying fhe frigger selecfion are required fo have af leasf one reconstiucfed vertex wifh af leasf 
five associafed fracks wifh px > 400 MeV, consisfenf wifh originafing from fhe beam collision region in 
fhe x-y plane. The average number of pp inferacfions per bunch crossing is approximafely 20, resulfing in 
several vertices reconstiucfed per evenf. If more fhan one vertex is found, fhe hard-scaffer primary vertex 
is faken fo be fhe one which has fhe largesf sum of fhe squared fransverse momenfa of ifs associafed fracks. 
For fhe evenf fopologies considered in fhis paper, fhis requiremenf leads fo a probabilify fo reconsfrucf 
and selecf fhe correcf hard-scaffer primary vertex larger fhan 99%. 

Preselecfed evenfs are required fo have exacfly one election or muon, as defined in secfion 3, fhaf mafches, 
wifhin AT? = 0.15, fhe lepfon candidafe reconstiucfed by fhe frigger. In addifion, af leasf four jefs are 
required, of which af leasf fwo musf be T^-fagged. 


Light-jet denotes a jet originating from the hadronisation of a light quark (u, d, s) or gluon. 
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5. Background and signal modelling 


After the event preselection, the main background is tf WbWb production, possibly in association 
with jets, denoted by tt+jets in the following. Single top quark production and production of a Vk boson 
in association with jets (Vk+jets) contribute to a lesser extent. Small contributions arise from multijet, 
Z+jets and diboson {WW, fkZ, ZZ) production, as well as from the associated production of a vector 
boson V {V = W, Z) or a Higgs boson and a tt pair {ttV and ttH). Signal and all backgrounds are 
estimated from simulation and normalised to their theoretical cross sections, with the exception of the 
multijet background, which is estimated with data-driven methods [39]. 

Simulated samples of tt events are generated with the NLO generator Powheg-Box 2.0 [40-43] using 
the CTIO [44] set of parton distribution functions (PDF). The nominal sample is interfaced to Pythia 
6.425 [45] for parton showering and hadronisation with the CTEQ6L1 PDF set and the Perugia2011C [46] 
set of optimised parameters for the underlying event (UE) description, referred to as the “UE tune”. An al¬ 
ternative sample, used to study the uncertainty related to the hadronisation model, is interfaced to Herwig 
v6.520 [47] with the CTEQ6E1 PDE set and Jimmy v4.31 [48] to simulate the UE. All samples are gener¬ 
ated assuming a top quark mass of 172.5 GeV and top quark decays exclusively through t Wb. The tt 
process is normalised to a cross section of 253];}g pb, computed using Top-i-e v2.0 [49] at next-to-next-to- 
leading order (NNEO) in QCD, including resummation of next-to-next-to-leading logarithmic (NNEE) 
soft gluon terms [50-54], and using the MSTW 2008 NNEO [55, 56] PDE set. Theoretical uncertainties 
result from variations of the factorisation and renormalisation scales, as well as from uncertainties on the 
PDE and as- The latter two represent the largest contribution to the overall theoretical uncertainty on 
the cross section and were calculated using the PDE4EHC prescription [57] with the MSTW 2008 68% 
CE NNEO, CTIO NNEO [44, 58] and NNPDE2.3 5f FFN [59] PDE sets. In the case where a non-zero 
BR(f ^ Hq) is assumed, an additional factor of [1 - BR(t —> Hq)]^ is applied to the sample normalisation. 
It is not possible to generate the ti —> WbHq signal with Powheg-Box, and a different event generator is 
used instead, as discussed below. 

The tt samples are generated inclusively, but events are categorised depending on the flavour content of 
additional particle jets not originating from the decay of the ti system.^ Details about this categorisation 
scheme can be found in ref. [21]. In this way, a distinction is made between ti+bb, ti+cc and tt-flight-jets 
events. The first two categories are generically referred to as tt-fHE events (with HE standing for “heavy 
flavour”), while the latter category also includes events with no additional jets. The modelling of tt-fHE in 
Powheg-Box- fPYTHiA is via the parton-shower evolution. To study uncertainties related to this simplified 
description, an alternative tt-fjets sample is generated with Madgraph5 1.5.11 [60] using the CTIO PDE 
set. It includes tree-level diagrams with up to three additional partons (including b- and c-quarks) and is 
interfaced to Pythia 6.425. 

Since the best possible modelling of the K-fjets background is a key aspect of this search, a correction is 
applied to simulated ti events in Powheg-Box-ePythia based on the ratio of the differential cross sections 
measured in data and simulation ai ^fs = 1 TeV as a function of top quark pj and ti system pj [61]. This 
correction significantly improves agreement between simulation and data at = 8 TeV in distributions 
such as the jet multiplicity and the px of decay products of the ti system [21], and is applied only to 
tt-flight-jets and ti + cc events. The modelling of the ti + bb background is improved by reweighting 
the Powheg-Box -fPYTHiA prediction to an NEO prediction of ti + bb with massive b quarks and including 


^ Particle jets are reconstructed by clustering stable particles excluding muons and neutrinos using the anti-fe, algorithm with a 
radius parameter R = 0.4. 
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parton showering [62], based on Sherpa+OpenLoops [63, 64] using the CTIO PDF set. Such treatment is 
not possible for the tt + cc background since a corresponding NLO prediction is not currently available. 
More details about the modelling of the tf+jets background can be found in ref. [21]. 

Samples of single-top-quark backgrounds corresponding to the t-channel, 5-channel, and Wt production 
mechanisms are generated with Powheg-Box 2.0 [65, 66] using the CTIO PDF set and interfaced to Pythia 
6.425 with the CTEQ6L1 PDF set in combination with the Perugia2011C UE tune. Overlaps between the 
tt and Wt final states are avoided using the “diagram removal” scheme [67]. The single-top-quark samples 
are normalised to the approximate NNEO theoretical cross sections [68-70], calculated using the MSTW 
2008 NNEO PDE set. 

Samples of W/Z-i-jets events are generated with up to five additional partons using the Alpgen v2.14 [71] 
EO generator with the CTEQ6E1 PDE set and interfaced to Pythia v6.426. To avoid double-counting of 
partonic configurations generated by both the matrix-element calculation and the parton shower, a parton- 
jet matching scheme (“MEM matching”) [72] is employed. The VF-tjets samples are generated separately 
for VF-rlight-jets, Wbb+]&t?,, fFcc-i-jets, and ITc-f-jets. The Z-i-jets samples are generated separately for 
Z-i-light-jets, ZfiD-i-jets, and Zcc-i-jets. Overlap between VQQ+]e.t?> {V = W,Z and Q = b,c) events 
generated from the matrix-element calculation and those generated from parton-shower evolution in the 
IF/Z-Elight-jets samples is avoided via an algorithm based on the angular separation between the extra 
heavy quarks: if AR{Q, Q) > 0.4, the matrix-element prediction is used, otherwise the parton-shower 
prediction is used. Both the fF-i-jets and Z-i-jets background contributions are normalised to their inclusive 
NNEO theoretical cross sections [73]. Eurther corrections are applied to W/Z+jets events in order to 
better describe data in the preselected sample. Normalisation factors for each of the IT-i-jets categories 
{Wbb+iets, ITcc-i-jets, ITc-f-jets and IT-i-light-jets) are derived for events with one lepton and at least four 
jets by simultaneously analysing six different event categories, defined by fhe b-ta.g mulfiplicify (0, 1 
and >2) and fhe sign of fhe lepfon charge [74]. The b-tag mulfiplicify provides informafion abouf fhe 
heavy-flavour composifion of fhe IT-i-jefs background, while fhe lepfon charge is used fo defermine fhe 
normalisafion of each componenf, exploifing fhe expecfed charge asymmefry for VF-i-jefs producfion in pp 
collisions as predicfed by Alpgen. In fhe case of Z-i-jels evenfs, a correcfion fo fhe heavy-flavour fraclion is 
derived fo reproduce fhe relative rales of Z-i-2-jels evenfs wilh zero and one fi-lagged jef observed in dala. 
In addifion, fhe Z boson pj specfrum is compared befween dala and fhe simulalion in Z-i-2-jels evenfs, 
and a reweighling funclion is derived in order fo improve fhe modelling. This reweighling funcfion is 
also applied fo fhe VF-tJefs simulaled sample and if was verified lhal Ihis correcfion further improves fhe 
agreemenf befween dafa and simulalion for IF-i-jels evenfs. In any case, IF/Z-i-jefs evenfs conslilule a very 
small background in fhis analysis affer final evenl selecfion. 

The VFVF/VFZ/ZZ-pjels samples are generated wilh up fo Ihree additional parlons using Alpgen v2.13 and 
fhe CTEQ6E1 PDF sel, inferfaced fo Herwig v6.520 and Jimmy v4.31 for parton showering, hadronisalion 
and UE modelling. The MEM parfon-jel malching scheme is used. The IFIF-i-jefs samples require al leasl 
one of fhe W bosons fo decay leplonically, while fhe IFZ/ZZ-i-jefs samples require one Z boson fo decay 
leplonically and olher boson decays inclusively. Additionally, IFZ-i-jefs samples requiring fhe W boson 
fo decay leplonically and fhe Z boson fo decay hadronically, are generaled wilh up fo Ihree additional 
parlons (including massive b- and c-quarks) using Sherpa v 1.4.1 and fhe CTIO PDF sef. All diboson 
samples are normalised fo fheir NEO Iheorelical cross seclions [75]. 

Samples of ttV evenfs, including ttWW, are generaled wilh up fo Iwo addilional parlons using Madgraph5 
1.3.28 wilh fhe CTEQ6E1 PDF sel, and inferfaced fo Pythia 6.425 wilh fhe AUET2B UE lune [76]. A 
sample of ttH evenfs is generated wilh fhe PowHel framework [77], which combines fhe Powheg-Box 


generator and NLO matrix elements obtained from the HELAC-Oneloop package [78]. The sample is gen¬ 
erated using the CTlOnlo PDF set [44]. Showering is performed with Pythia 8.1 [79] using the CTEQ6L1 
PDF set and the AU2 UE tune [76, 80]. Inclusive decays of the Higgs boson are assumed in the generation 
of the ttH sample. The tW samples are normalised to the NEO cross-section predictions [81]. The ttH 
sample is normalised using the NEO cross section [82-84] and the Higgs decay branching ratios [85-88] 
collected in ref. [89]. 

The multijet background contributes to the selected data sample via several production and misrecon- 
struction mechanisms. In the electron channel, it consists of non-prompt electrons (from semileptonic 
b- or c-hadron decays) as well as misidentified photons (e.g. from a conversion of a photon into an e'^e~ 
pair) or jets with a high fraction of their energy deposited in the EM calorimeter. In the muon chan¬ 
nel, the multijet background is predominantly from non-prompt muons. Its normalisation and shape 
are estimated directly from data by using the “matrix method” technique [39], which exploits differ¬ 
ences in lepton-identification-related properties between prompt and isolated leptons and leptons that are 
either non-isolated or result from the misidentification of photons or jets. Further details can be found in 
ref. [22]. 

The tt —> WbHq signal process is modelled using the Protos v2.2 [90, 91] EO generator with the 
CTEQ6E1 PDF set, and interfaced to Pythia 6.426 and the Perugia2011C UE tune. Two separate samples 
are generated corresponding to tt WbHc and tt WbHu, with the W boson forced to decay lepton- 
ically, W ^ bv {€ = e,ju,T), The top quark and Higgs boson masses are set to 172.5 GeV and 125 GeV, 
respectively. The Higgs boson is allowed to decay to all SM particles with branching ratios as given in 
ref. [89]. The signal sample is normalised to the same NNEO cross section as used for the ti WbWb 
sample, and the corresponding branching ratios: cr{tt W{-^ £v)bHq) = 2BR(t ^ Hq)[l - BR(t ^ 
//^)]BR(IT ^ bv)crth with BR(VT ^ £v) = 0.324 and BR(t ^ Hq) depending on the branching ratio 
being tested. Typically a reference branching ratio of BR(t —> Hq) = 1% is used. The case of both top 
quarks decaying into Hq is neglected in the analysis given existing upper limits on BR(f —> Hq) (see sec¬ 
tion 1). In order to improve the modelling of the signal kinematics, a two-step reweighting procedure is 
applied: the first step is designed to correct the spectrum of top quark pj and tt system pj to match that of 
the uncorrected tt —> WbWb Powheg-Box-i-Pythia sample; the second step involves the same correction 
to the top quark px and tt system px applied to the tf-i-jets background (see discussion above). 

Finally, all generated samples are processed through a simulation [92] of the detector geometry and 
response using Geant4 [93]. Additional minimum-bias pp interactions are simulated with the Pythia 
8.1 generator with the MSTW 2008 EO PDF set and the A2 UE tune [94]. They are overlaid on the 
simulated signal and background events according to the luminosity profile of fhe recorded dafa. The 
confribufions from fhese pileup inferacfions are modelled bofh wifhin fhe same bunch crossing as fhe hard- 
scalfering process and in neighbouring bunch crossings. All simulafed samples are processed fhrough fhe 
same reconsfrucfion soffware as fhe dafa. Simulated evenfs are correcfed so fhaf fhe objecf idenfificalion 
efficiencies, energy scales, and energy resolufions mafch fhose defermined from dafa confrol samples. 


6. Analysis strategy 

This section presenfs an overview of fhe analysis sfrafegy followed by fhe tt —> WbHq, H ^ bb search. 
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6.1. Event categorisation 

Given the focus on the W ^ (v and H ^ bb decay modes, the tt —> WbHq signal is expected to have 
typically four jets, of which three or four are Zj-tagged. The latter case corresponds to the tt WbHc 
signal where the charm quark, as well as the three Zj-quark jets, are Z>-tagged. Additional jets can also be 
present because of initial- or final-state radiation. In order to optimise the sensitivity of the search, the 
selected events are categorised into different channels depending on the number of jets (4, 5 and >6) and 
on the number of f^-tagged jets (2, 3 and >4). Therefore, the total number of analysis channels considered 
in this search is nine: (4 j, 2 b), (4 j, 3 b), (4 j, 4 b), (5 j, 2 b), (5 j, 3 b), (5 j, >4 b), (>6 j, 2 b), (>6 j, 3 b), 
and (>6 j, >4 b), where (n j, m b) indicates n selected jets and m fj-tagged jets. 

The overall rate and composition of the tf+jets background strongly depends on the jet and b-ta.g multi¬ 
plicities, as illustrated in figure 1. The tt-i-lighf-jefs background is dominanf in evenfs wifh exacfly fwo or 
fhree f?-fagged jefs, wifh fhe fwo f?-quarks from fhe fop quark decays being tagged in both cases, and a 
charm quark from the hadronic W boson decay also being tagged in the latter case. Contributions from 
tt+cc and tt+bb become significant as the jet and b-tag multiplicities increase, with the tt+bb background 
being dominant for events with >6 jets and >4 b-tags. 

In the channels with four or five jefs and fhree or af leasf four f7-fags, which dominafe fhe sensifivify of 
fhis search, selecfed signal evenfs have a H ^ bb decay in more fhan 95% of fhe evenfs. The channels 
mosf sensifive fo fhe tt —> WbHu and tt —> WbHc signals are (4 j, 3 b) and (4 j, 4 b) respecfively. 
Because of fhe beffer signal-fo-background rafio in fhe (4 j, 4 b) channel, fhis analysis is expecfed fo 
have beffer sensifivify for tt WbHc fhan for ti WbHu signal. The resf of fhe channels have 
significanlly lower signal-fo-background rafios, buf fhey are useful for calibrating fhe tf-i-jels background 
prediction and consfraining fhe related sysfemafic uncerfainfies (see secfion 7) fhrough a likelihood fif fo 
dafa (see secfion 8). This sfrafegy was firsf used in fhe ATLAS search for ttH associafed production, wifh 
H —> bb [21], and is adopfed in fhis analysis. A fable summarising fhe observed and expecfed yields 
before fhe til fo dafa in each of fhe analysis channels can be found in appendix A. 


6.2. Discrimination of signal from background 

Afler even! cafegorisafion, fhe signal-fo-background rafio is very low even in fhe mosf sensifive analysis 
channels, and a suifable discriminating variable befween signal and background needs fo be consfrucfed 
in order fo improve fhe sensifivify of fhe search. A powerful discriminanf befween signal and background 
can be defined as: 

D{x) - —^-jT-, (2) 

where P®‘®(x) and P'’'^®(x) represenf fhe probabilify densify funcfions (pdf) of a given evenf under fhe sig¬ 
nal hypofhesis {tt WbHq) and under fhe background hypofhesis {tt WbWb) respectively. Bofh pdfs 
are functions of x, representing fhe four-momenfum vecfors of all final-sfafe particles af fhe reconsfrucfion 
level: fhe lepfon {{), fhe neufrino (v; reconsfrucfed as discussed below), and fhe Ajets selecfed jefs in a 
given analysis channel. 

Since bofh signal and background resulf from fhe tt decay, fhere are few experimenfal handles available 
fo discriminate befween fhem. The mosf prominenf feafures are fhe differenf resonances presenf in fhe 
decay (i.e. fhe Higgs boson in the case of tt —> WbHq and a hadronically decaying W boson in the case 
of ti WbWb), and the different flavour content of the jets forming those resonances. This is the main 
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Figure 1: Comparison between the data and background prediction for the yields in each of the analysis channels 
considered before the fit to data (pre-fit). Backgrounds are normalised to their nominal cross sections discussed in 
section 5. The expected ff —> WbHc and tt —> WbHu signals (dashed histograms) are shown separately normalised 
to BR(f —> Hq) - 1%. The tt WbWb background is normalised to the SM prediction. The small contributions 
from VT/Z-i-jets, single top, diboson and multijet backgrounds are combined into a single background source referred 
to as “Non-ff”. The bottom panel displays the ratio of data to the SM background (“Bkg”) prediction. The hashed 
area represents the total uncertainty on the background. 


information exploited in the construction of and P’’’^§(x) in this analysis, so that x is extended to 

include not only the four-momenta of jets pjet, but also the value of their multivariate ^-tagging discrimi¬ 
nant Wjet, i.e., X = {p(, py, (pjet,, (* = 1, • • • , A^jets)- There is also some angular information from the 
different spins of the daughter resonances (Higgs and W boson) that could be exploited, but it is expected 
to be subleading in importance and is neglected in this analysis. 

The calculation of P'^‘“(x) and P^'^“(x) is discussed in detail in sections 6.2.1 and 6.2.2 respectively. In 
the following, denotes the A-quark jet from the semileptonic top quark decay, qn and bh denote the 
light-quark jet {qu - u or c) and A-quark jet from the hadronic top quark decay in background and signal 
events respectively, q\ and <72 denote the up-type-quark jet (u or c) and down-type-quark jet {d or s) from 
the W boson decay respectively, and bi and A 2 denote the two A-quark jets from the Higgs boson decay. 
The level of separation achieved between signal and background with the resulting discriminant D is 
illustrated in section 6.2.3. 


6.2.1. Signal probability 

The construction of P®‘®(x) will now be described step by step to illustrate the method. If the partonic 
origin of each jet were known [see figure 2 (a)], P'^‘S(x) would be defined in fhis analysis as fhe producf 
of fhe normalised pdfs for each of fhe reconsfrucfed invarianf masses in fhe evenf: fhe semilepfonic fop 
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quark mass the hadronic top quark mass {Mb^h^qj^) and the Higgs boson mass Since 

Mbfb 2 qh Mbib 2 correlated, their difference in quadrature, Xb^b 2 qh = ^bib 2 qh instead 

of Mb^b 2 qh- Therefore the expression for just making use of the above kinematic information, denoted 



P^f„(x) = P^HMme)P"KXb,b2q,)P"KMb,b2)- 


(3) 


The distributions of these invariant masses are obtained from simulated signal events using the recon¬ 
structed lepton and/or jets corresponding to the correct parton-jet assignment, determined by matching 
a given quark (before final-state radiation) to the closest jet with AT? < 0.3. The corresponding pdfs 
are constructed as unit-normalised one-dimensional histograms. To compute M^ybf, the neutrino four- 
momentum is needed, which is reconstructed as follows. Initially, the x and y components of the neutrino 
momentum, y and py^y, are identified wifh fhose of fhe reconsfrucfed vector. The z component of 
the neutrino momentum, p^ y, is inferred by solving = (p^ -i- Pv)^, with = 80.4 GeV being the 
W boson mass. If two real solutions (“2sol”) exist, they are sorted according to their absolute value of 
|pj,y| i.e., |pz,yi| < \Pz,v 2 \- It is found that in 62% of the cases p^^yi is closer than p^^yi to the generator-level 
neutrino Pz,v In this case, two different pdfs are constructed, one for each solution, and P^^^^^{Mbybf) is 
defined as fhe average of fhe two pdfs weighted by their fractions (0.62 for p^ yi and 0.38 for Pz,v 2 )- If 
no real solution (“nosol”) exists, which happens in about 30% of the cases, the px,v and py^y components 
are scaled by a common factor until the discriminant of the quadratic equation is exactly zero, yielding 
only one solution for p^ y. This solution for p^ y is used to compute Mtyb^, from which the correspond¬ 
ing P^^l^aii^tvbt) is constructed. In the calculation of T’y®(x) from equation (3), P^^^{M(ybi) is identified 
wifh P^ 2 ^f^\{Mtybf) or wifh P^^^^^^{Mevbf), depending on how many neutrino solutions can be found for the 
event. 

In practice, the partonic origin of the jets is not known, so it is necessary to evaluate P®‘®(x) by averaging 
over the Np possible parton-jet assignments, which dilutes the kinematic information. At this point b- 
tagging information can be used to suppress the impact from parton-jet assignments that are inconsistent 
with the correct parton flavours as follows: 



(4) 


where T’^;^(x) is given by equafion (3) and T’^(®g(x) is defined as: 

^btagW = ■Pfo(jeti)/’fe(jet2)/’fe(jet3)^<?/.(jet4)> 


(5) 


wifh jef,- {i - 1,... ,4) represenfing fhe parfon-jef assignment being evaluated, and Py(jet,) denoting the 
probability that jet i, characterised by its four-momentum pjet. and T^-tagging weight value Wjet., originates 
from a parton with flavour / {b, c, or /; I for light parton). The calibration of the /^-tagging algorithm 
is performed for fixed thresholds on the multivariate /7-tagging discriminant variable, corresponding to 
different average /7-tagging efficiencies in tt events of 60%, 70%, and 80%, also referred to as “operating 
points” (OP). The corresponding thresholds are denoted by with OP = 60%, 70%, or 80%. Parame- 
terisations of the /7-tagging efficiencies for different jet flavours as functions of jet pj and p are available 
for each of these operating points, e^^{pT:, p), which can be used to compute Pf as follows: if the jet b- 

tagging weight falls between the thresholds for operating points OPi and OP 2 , < wjet < then 
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Pf = 6°^' - ; alternatively, if the jet Z>-tagging weight is below (above) the threshold corresponding 

to the 80% (60%) operating point, then Pf - I - {Pf - 
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Figure 2: Representative Feynman diagrams illustrating the partonic configurations and parton-jet assignments 
considered in the construction of (a) the signal probability and (b) and (c) the background probability used in the 
definition of the final discriminant (see text for details). 


6.2.2. Background probability 

The calculation of P^^^ follows a similar approach to that discussed in section 6.2.1, although it is slightly 
more complicated to account for the varying fraction and different kinematic features of the ff+light-jets, 
tt + cc and tt + bb backgrounds as a function of the analysis channel. This is particularly relevant in the 
(4 j, 3 b) and (4 j, 4 b) channels, which dominate the sensitivity of the search. While ff+light-jets events 
often have both jets from the hadronic W boson decay among the four selected jets [see figure 2(b)], this 
is seldom the case for tt + bb and tt + cc events, especially in the (4 j, 4 b) channel. In this case the 
four Zj-tagged jets typically originate from the two Z^-quarks from the top quark decays, the charm quark 
from the W boson decay, and an extra heavy-flavour quark {b or c) produced in association with the tt 
system, while the jet associated with the down-type quark from the W boson decay is not reconstructed 
[see figure 2(c)]. 

To account for this, the following kinematic variables are considered: Xqjth Mqj, with 

Xqjhi, = Xlq^jhi^ 0 Mq^j, wcrc j dcnotcs an extra quark-jet which can either originate from the W bo- 
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son decay (^ 2 ) or from an extra heavy-quark {b or c) produced in association with the tt system. For each 
of these possibilities, occurring in a fraction fj of the cases, corresponding pdfs are constructed. As a 
generalisation of equation (4), the expression for becomes: 


pbkg(^) 


z 

k=\ 


I. 

je{b,c,q2] 


f .pbkg,; 
J } btag 


(x^)F 


bkg,i 

kin 


(x') 


'V 

i: 

k=\ 


z 

je[b,c,q 2 \ 


f .pbkg,;' 
J j btag 


(X^) 


( 6 ) 


with 

and 

^Sg'W = ^&(ieti)P,, aet2)Pyaet3)Pz,aet4). (8) 

where P/^(jet,) are computed as discussed in section 6.2.1. In the above expression, Pj = Pi for j - < 72 , 
the down-type quark in the W boson decay, and Pq^ = fcPc + (1 - fc)Pu where fc is the fraction of events 
where the up-type quark from the W boson decay assigned to the jet is a charm quark. This fraction is 
different in each analysis channel, primarily depending on the Z^-tag multiplicity requirements. It varies 
from ~ 50% for events in the (4 j, 2 b) channel to ~ 90% for events in the (4 j, 4 b) channel. 


6.2.3. Final discriminant 

The final discriminant D is computed for each event as given in equation (2), using the definitions for P^^^ 
and given in equations (4) and (6), respectively. Since this analysis has higher expected sensitivity 
to a ft —> WbHc signal than to att ^ WbHu signal and, in order to allow probing of the BR(f —> PIu) 
versus BR(t —> He) plane, the discriminant optimised for tt —> WbHc is used for both the He and Hu 
decay modes. It was verified that using the tt —> WbHc discriminant for the tt —> WbHu search does 
not result in a significant sensitivity loss. Figure 3 compares the shape of the D distribution between the 
tt WbHc and tt WbHu signals and the tt WbWb background in each of the channels considered 
in this analysis. 


7. Systematic uncertainties 

Several sources of systematic uncertainty are considered that can affect the normalisation of signal and 
background and/or the shape of their corresponding final discriminant distributions. Each source of sys¬ 
tematic uncertainty is considered to be uncorrelated with the other sources. Correlations of a given sys¬ 
tematic uncertainty are maintained across processes and channels. Table 1 presents a list of all systematic 
uncertainties considered in the analysis and indicates whether they are taken to be normalisation-only, or 
to affect both shape and normalisation. 

The leading sources of systematic uncertainty vary depending on the analysis channel considered, but 
they typically originate from ff-i-jets modelling (including tf-i-HF) and fj-tagging. For example, the total 
systematic uncertainty in the background normalisation in the (4 j, 4 b) channel, which dominates the 
sensitivity in the case of the tt —> WbHc search, is approximately 20%, with the largest contributions ori¬ 
ginating from tt-i-HF normalisation, / 7 -tagging efficiency, c-tagging efficiency, light-jet tagging efficiency 
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(c) 





(e) 
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(g) 


(h) 


(i) 


Figure 3; Comparison of the shape of the D discriminant distribution between the tt —> WbHc (red dashed) and 
tt WbHu (blue dotted) signals, and the tt WbWb background (black solid) in each of the channels considered 
in the analysis: (a) (4 j, 2 b), (b) (4 j, 3b), (c) (4 j, 4 b), (d) (5 j, 2 b), (e) (5 j, 3 b), (f) (5 j, >4 b), (g) (>6 j, 2 b), (h) 
(>6j, 3 b), and(i) (>6j, >4 b). 
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Systematic uncertainty 

Type 

Components 

Luminosity 

N 

1 

Reconstructed Objects 



Electron 

SN 

5 

Muon 

SN 

6 

Jet reconstruction 

SN 

1 

Jet vertex fraction 

SN 

1 

Jet energy scale 

SN 

22 

Jet energy resolution 

SN 

1 

Missing transverse momentum 

SN 

2 

^-tagging efficiency 

SN 

6 

c-tagging efficiency 

SN 

4 

Light-jet tagging efficiency 

SN 

12 

High-;rT tagging 

SN 

1 

Background Model 



tt cross section 

N 

1 

tt modelling: pj reweighting 

SN 

9 

tt modelling: parton shower 

SN 

3 

tt+HF: normalisation 

N 

2 

tt+cc\ Pi reweighting 

SN 

2 

tt+cc: generator 

SN 

4 

tt+bb: NLO shape 

SN 

8 

VE-i-jets normalisation 

N 

3 

W Pi reweighting 

SN 

1 

Z-i-jets normalisation 

N 

3 

Z Pi reweighting 

SN 

1 

Single top normalisation 

N 

3 

Single top model 

SN 

1 

Diboson normalisation 

N 

3 

ttV cross section 

N 

1 

ttV model 

SN 

1 

ttH cross section 

N 

1 

ttH model 

SN 

2 

Multijet normalisation 

N 

4 

Signal Model 



tt cross section 

N 

1 

Higgs boson branching ratios 

N 

3 

tt modelling: pi reweighting 

SN 

9 

tt modelling: pi reweighting non-closure 

N 

1 

tt modelling: parton shower 

N 

1 


Table 1: List of systematic uncertainties considered. An “N” means that the uncertainty is taken as affecting only 
the normalisation for all relevant processes and channels, whereas “SN” means that the uncertainty is taken on both 
shape and normalisation. Some of the systematic uncertainties are split into several components for a more accurate 
treatment. 
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and tt cross section. However, as shown in section 9, the fit to data in the nine analysis channels allows 
the overall background uncertainty to be reduced significantly, to approximately 4.4%. The reduced un¬ 
certainty results from the significant constraints provided by the data on some systematic uncertainties, 
as well as the anti-correlations among sources of systematic uncertainty resulting from the fit to the data. 
The total systematic uncertainty on the tt WbHc signal normalisation in the (4 j, 4 b) channel is 
approximately 17%, with similar contributions from uncertainties related to f^-tagging and overall signal 
modelling. After the fit, this uncertainty is reduced to 7.8%. Table 2 presents a summary of the systematic 
uncertainties for the tt —> WbHc search and their impact on the normalisation of the signal and the main 
backgrounds in the (4 j, 4 b) channel. 

The following sections describe each of the systematic uncertainties considered in the analyses. 



WbHc 

Pre-fil 

t7±LJ tt + cc 

tt + bb 

WbHc 

Post-fil 

t7±LJ tt + cc 

tt + bb 

Luminosity 

±2.8 

±2.8 

±2.8 

±2.8 

±2.6 

±2.6 

±2.6 

±2.6 

Lepton efficiencies 

±1.5 

±1.5 

±1.5 

±1.5 

±1.5 

±1.5 

±1.5 

±1.5 

Jef energy scale 

±3.3 

±2.9 

±2.3 

±5.8 

±1.4 

±1.2 

±1.8 

±4.1 

Jef efficiencies 

±1.2 

- 

±1.9 

±1.7 

±0.9 

- 

±1.4 

±1.2 

Jef energy resolufion 

- 

±1.2 

±2.8 

±2.9 

- 

- 

±1.0 

±1.1 

f>-fagging eff. 

±7.9 

±5.5 

±5.2 

±10 

±5.7 

±3.9 

±3.7 

±6.6 

c-fagging eff. 

±7.0 

±6.6 

±13 

±3.5 

±6.3 

±6.0 

±11 

±3.2 

Lighf-jef fagging eff. 

±0.8 

±18 

±3.2 

±1.5 

±0.6 

±13 

±2.3 

±1.1 

tt: reweighfing 

±5.9 

±2.7 

±4.2 

- 

±3.8 

±1.9 

±2.3 

- 

tt: parfon shower 

±5.4 

±4.8 

±10 

±4.9 

±1.7 

±1.5 

±6.5 

±3.1 

?7-i-HF: normalisafion 

- 

- 

±50 

±50 

- 

- 

±32 

±16 

tf-i-HF: modelling 

- 

- 

- 

±7.7 

- 

- 

- 

±7.4 

Signal modelling 

±6.9 

- 

- 

- 

±6.9 

- 

- 

- 

Theor. cross secfions 

±6.2 

±6.2 

±6.2 

±6.2 

±3.9 

±3.9 

±3.9 

±3.9 

Tofal 

±17 

±22 

±54 

±53 

±7.8 

±14 

±28 

±15 


Table 2: tt —> WbHc,H —> bb search: summary of the systematic uncertainties considered in the (4 j, 4 b) channel 
and their impact (in %) on the normalisation of the signal and the main backgrounds, before and after the fit to 
data. The tt —> WbHc signal and the ff-rlight-jets background are denoted by “WbHc"" and “tT+LJ"" respectively. 
Only sources of systematic uncertainty resulting in a normalisation change of at least 0.5% are displayed. The total 
post-fit uncertainty can differ from the sum in quadrature of individual sources due to the anti-correlations between 
them resulting from the fit to the data. 


7.1. Luminosity 

The uncertainty on the integrated luminosity is 2.8%, affecting the overall normalisation of all processes 
estimated from the simulation. It is estimated from a calibration of the luminosity scale derived from 
beam-separation scans performed in November 2012, following the same methodology as that detailed in 
ref. [38]. 
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7.2. Reconstructed objects 


Uncertainties associated with leptons arise from the reconstruction, identification and trigger. These effi¬ 
ciencies are measured using tag-and-probe techniques on Z —> {{ - e,ii) data and simulated samples. 

The small differences found are correcfed for in fhe simulation. Negligible sources of uncerfainfy origin- 
ale from fhe correclions applied fo adjusf fhe lepfon momenlum scale and resolufion in fhe simulation 
lo mafch Ihose in dafa. The combined effecl of all fhese uncerlainfies resulls in an overall normalisation 
uncertainly on fhe signal and background of approximalely 1.5%. 

Uncerlainfies associated wilh jels arise from fhe efficiency of jel reconslrucfion and idenlificalion based 
on fhe JVF variable, as well as fhe jel energy scale and resolution. The largesl conlribufion resulls from 
fhe jel energy scale, whose uncerlainly dependence on jel px and rj is splil into 22 uncorrelaled sources 
lhal are Irealed independenlly in Ihe analysis. If affecls Ihe normalisation of signal and backgrounds by 
approximately 3^% in Ihe mosl sensitive search channels, (4 j, 3 b) and (4 j, 4 b), and up to 12% in Ihe 
channels wilh >6 jels. 

Uncerlainlies associated wilh energy scales and resolutions of leptons and jels are propagated to F™®®. 
Additional uncertainties originating from Ihe modelling of Ihe underlying evenl, in particular ils impacl 
on Ihe px scale and resolution of uncluslered energy, are negligibly small. 

The leading uncerlainfies associated wilh reconslrucled objecls in Ihis analysis originate from Ihe mo¬ 
delling of Ihe b-, C-, and lighl-jel-lagging efficiencies in Ihe simulation, which is corrected to malch Ihe 
efficiencies measured in dala conlrol samples [36, 95] Ihrough dedicated scale factors. Uncertainties 
on Ihese factors include a lolal of six independenl sources affecting l?-jels and four independenl sources 
affecting c-jels. Each of Ihese uncertainties has a differenl jel-px dependence. Twelve sources of uncer¬ 
lainly affecting lighl jels are considered, which depend on jel px and p. The above sources of systematic 
uncertainly are laken as uncorrelaled belween l^-jels, c-jels, and lighl-jels. They have Iheir largesl impacl 
in Ihe (4 j, 4 b) channel, resulting in 10%, 13%, and 18% normalisation uncertainties on Ihe ti + bb, 
tt + cc, and rt+lighl-jels background associated wilh Ihe uncerlainfies on Ihe b-, c-, and lighl-jel-lagging 
scale factors, respectively. An additional uncerlainly is included due to Ihe exlrapolalion of Ihese scale 
factors to jels wilh px beyond Ihe kinematic reach of Ihe dala calibration samples used (px > 300 GeV 
for b- and c-jels, and px > 750 GeV for lighl-jels), laken to be correlated among Ihe Ihree jel flavours. 
This uncerlainly has a very small impacl in Ihis analysis (e.g. < 0.2% on Ihe signal and background 
normalisations in Ihe (4 j, 4 b) channel). 


7.3. Background modelling 

A number of sources of systematic uncerlainly affecting Ihe modelling of tf+jels are considered. A 
brief summary is provided below, wilh furlher delails available in ref. [21], as Ihe uncerlainly Irealmenl is 
identical. An uncerlainly of -i-6.1%/-6.4% is assumed for Ihe inclusive tt production cross section [49], in¬ 
cluding conlribulions from varying Ihe factorisation and renormalisalion scales, and uncerlainfies arising 
from Ihe PDF, as, and Ihe top quark mass. Uncerlainfies associated wilh Ihe reweighling procedure ap¬ 
plied to ft-flighl-jels and tt+cc processes include Ihe nine leading sources of uncerlainly in Ihe differential 
cross-section measuremenl at ^fs = 1 TeV [61]. Additional uncertainties assigned to Ihe modelling of Ihe 
tt + cc background include a 50% normalisation uncerlainly, Ihe full differences belween applying and nol 
applying Ihe reweighlings of Ihe top quark and tt px speclra, as well as smaller uncerlainfies associated 
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with the choice of LO generator. Uncertainties affecting the modelling of tt + bb production include a nor¬ 
malisation uncertainty of 50% (taken to be uncorrelated with the same uncertainty assigned to the tt + cc 
background) and shape uncertainties (including inter-category migration effects) associated with the NLO 
prediction from Sherpa-i-OpenLoops, which is used for reweighting of the default Powheg-Box tt + bb pre¬ 
diction. These include three different scale variations, a different shower-recoil model scheme, and two 
alternative PDF sets (MSTW and NNPDF). Additional uncertainties are assessed for the contributions to 
the tt+bb background originating from multiple parton interactions or final-state radiation from top decay 
products, which are not part of the NLO prediction. Finally, an uncertainty due to the choice of parton 
shower and hadronisation model is derived by comparing events produced by Powheg-Box interfaced to 
Pythia or Herwig. This uncertainty is taken to be uncorrelated between the tf-i-light-jets, tt + cc and tt + bb 
processes. 

Uncertainties affecting the modelling of the fP-i-jets background include a 7% normalisation uncertainty 
for events with >4 jets coming from the data-driven normalisation procedure. The corresponding normal¬ 
isation uncertainty for Z-i-jets is 5% for events with >2 jets. In addition, a 24% normalisation uncertainty 
is added in quadrature for each additional inclusive jet-multiplicity bin beyond the one where the back¬ 
ground is normalised, based on a comparison among different algorithms for merging LO matrix elements 
and parton shower simulations [96]. For example, IT-i-jets events with exactly 4 jets, exactly 5 jets and 
> 6 jets are assigned normalisation uncertainties of 7%, 7% 0 24% = 25% and 7% 0 24% 0 24% = 35%. 
Finally, the full size of the W and Z boson correction, after symmetrisation, is taken as a systematic 
uncertainty. The above uncertainties are taken as uncorrelated between IT-i-jets and Z-i-jets. 

Uncertainties affecting the modelling of the single-top-quark background include a -p 5%/-4% uncer¬ 
tainty on the total cross section, which is estimated as a weighted average of the theoretical uncertainties 
on t-, Wt- and ^-channel production [68-70]. Similarly to the case of IT/Z-pjets, an additional 24% nor¬ 
malisation uncertainty is added in quadrature for each additional inclusive jet-multiplicity bin above >3 
jets. An additional systematic uncertainty on ITf-channel production concerning the separation between 
tt and Wt at NLO [97] is assessed by comparing the nominal sample, which uses the so-called “diagram 
subtraction” scheme, with an alternative sample using the “diagram removal” scheme. 

Uncertainties on the diboson background normalisation include 5% from the NLO theoretical cross sec¬ 
tions [75] and additional 24% normalisation uncertainties added in quadrature for each additional in¬ 
clusive jet-multiplicity bin above >2 jets. Uncertainties on the ttV and tJH normalisations are 15% and 
-i-9%/-12% respectively, from the uncertainties on their respective NLO theoretical cross sections [81, 89, 
98]. Additional small uncertainties arising from scale variations, which change the amount of initial-state 
radiation and thus the event kinematics, are also considered. 

Uncertainties on the data-driven multijet background estimate receive contributions from the limited 
sample size in data, particularly at high jet and b-t&g multiplicities, as well as from the uncertainty on 
the rate of fake leptons, estimated in different control regions (e.g. selected with a requirement on either 
the maximum F"™*® or m^). A combined normalisation uncertainty of 50% due to all these effects is 
assigned, which is taken as correlated across jet and b-t&g multiplicity bins, but uncorrelated between 
electron and muon channels. No explicit shape uncertainty is assigned since the large statistical uncer¬ 
tainties associated with the multijet background prediction, which are uncorrelated between bins in the 
final discriminant distribution, effectively cover all possible shape uncertainties. 
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7.4. Signal modelling 


Several normalisation and shape uncertainties are taken into account for the tt —> WbHq signal. The 
uncertainty on the tt cross section (see above) also applies to the tt WbHq signal and is taken to be 
the same as, and fully correlated with, that assigned to the tt WbWb background. Uncertainties on the 
H ^ bb branching ratio are taken into account following the recommendation in ref. [89]: ±1.1% (Aas); 
±1.4% (Anih) and ±0.8% (theory). Additional modelling uncertainties originate from non-closure of the 
reweighting procedure applied to correct the distributions of top quark px and ti system px from Protos 
to match those from the uncorrected Powheg-Box±Pythia simulation, and the uncertainties associated 
with the further reweighting in the same variables to match the differential cross-section measurements 
at = 7 TeV, taken to be fully correlated with those assigned to the tt±light-jets background. Finally, 
an uncertainty from the choice of parton shower and hadronisation model is estimated by comparing the 
tt WbWb yields between Powheg-Box±Pythia and Powheg-Box±Herwig in the channels with two 
(7-tags, which are enriched in tf±light-jets, and thus taken to be representative of what would be the signal 
acceptance uncertainty due to differences in exfra jef radiafion and ( 7 -quark fragmenfafion befween fhe 
fwo parfon shower/hadronisafion models. 


8. Statistical analysis 


The disfribufions of fhe final discriminanfs in each of fhe analysis channels considered are combined fo 
fesf for fhe presence of a signal. The sfafisfical analysis is based on a binned likelihood funcfion L(p, 9) 
consfrucfed as a producf of Poisson probabilify ferms over all bins considered in fhe search. In fhe case 
of several searches being combined, fhe producf of Poisson probabilify ferms is exfended over all bins 
considered in all searches. The function L(p, 6) depends on fhe signal-sfrengfh paramefer p, defined as 
a mulfiplicafive facfor fo fhe yield for tt WbHq signal evenfs normalised fo a reference branching 
rafio BRref(t ^ Hq) = 1%, and 9, a sef of nuisance paramefers fhaf encode fhe effecf of sysfemafic 
uncerfainfies on fhe signal and background expecfafions and are implemenfed in fhe likelihood function 
as Gaussian or log-normal priors wifh fheir widfh parameters corresponding fo fhe size of fhe respecfive 
uncerfainfies. The relafionship befween p and fhe corresponding BR(t ^ Hq) is: 

^ BR(t ^ Hq)[\ - BR(t ^ Hq)] 

^ BR,ef(t ^ Hq){\ - BRref(t ^ Hq)]' ^ ’ 

For a given p value, fhe SM tt —> WbWb background confribufion is scaled accordingly in order fo 
preserve fhe inclusive tt cross secfion. The corresponding mulfiplicafive facfor would be [1 - BR(t ^ 
Hq)]^, wifh BR(t ^ Hq) being a funcfion of p as can be derived from equafion (9): 


BR(t 


Hq) = 


1 - Vl - 4BR,ef(f 


Hq){l - BRref(t 
2 


Hq))ti 


( 10 ) 


Therefore, fhe fofal number of signal and background evenfs in a given bin depends on p and 9. The 
besf-fif BR(t ^ Hq) is obfained by performing a binned likelihood fif fo fhe dafa under fhe signal- 
plus-background hypofhesis, i.e. maximising fhe likelihood funcfion L(p, 9) over p and 9. The nuisance 
paramefers 9 allow variafions of fhe expecfafions for signal and background according fo fhe corres¬ 
ponding sysfemafic uncerfainfies, and fheir tiffed values correspond fo fhe deviations from fhe nominal 
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expectations that globally provide the best fit to the data. This procedure allows a reduction of the im¬ 
pact of systematic uncertainties on the search sensitivity by taking advantage of the highly populated 
background-dominated channels included in the likelihood fit. 

The test statistic is defined as the profile likelihood ratio: = -21n(L(//, 0)), where ji and 0 

are the values of the parameters that maximise the likelihood function (with the constraint 0 < fi < ju), 
and are the values of the nuisance parameters that maximise the likelihood function for a given value 
of fi. Statistical uncertainties in each bin of the discriminant distributions are also taken into account 
via dedicated nuisance parameters in the fit. The test statistic is implemented in the RooFix pack¬ 
age [99, 100] and is used to measure the compatibility of the observed data with the background-only 
hypothesis by setting = 0 in the profile likelihood ratio: = -21n(L(0, The p-value 

(referred to as po) representing the compatibility of the data with the background-only hypothesis is es¬ 
timated by integrating the distribution of qo from background-only pseudo-experiments, approximated 
using the asymptotic formulae given in ref. [101], above the observed value of q^. The observed po"value 
is checked for each explored signal scenario. In the absence of any significant excess above the back¬ 
ground expectation, upper limits on p, and thus on BR(t ^ Hq) via equation (10), are derived by using 
q^ in the CLg method [102, 103]. Values of BR(t ^ Hq) yielding CLs<0.05, where CLg is computed 
using the asymptotic approximation [101], are excluded at >95% CL. 


9. Results 

This section presents the results obtained from the individual searches for tt —> WbHq, as well as their 
combination. 


9.1. H ^bb 

Following the statistical analysis discussed in section 8, a binned likelihood fit under the signal-plus- 
background hypothesis is performed on the distributions of the final discriminant in the nine analysis 
channels considered. Figures 4-6 show a comparison of the data and prediction in the final discriminant 
in each of the analysis channels, both pre- and post-fit to data, in the case of the tt WbHc search. 
The post-fit yields can be found in appendix A. The best-fit branching ratio obtained is BR(t ^ He) = 
[0.17 ± 0.12 (stat.) ± 0.17 (syst.)]%, assuming that BR(t ^ Hu) = 0. A similar fit is performed for 
the tt WbHu search, yielding BR(t —> Hu) = [-0.07 ± 0.17 (stat.) ± 0.28 (syst.)]%, assuming that 
BR(f —> He) = 0. The different measured values for the two branching ratios is the result of the different 
sensitivities of the tt WbHc and tt WbHu searches, as discussed in section 6.1. 

The large number of events in the analysis channels considered, together with their different background 
compositions, allows the fit to place constraints on the combined effect of several sources of systematic 
uncertainty. As a result, an improved background prediction is obtained with significantly reduced uncer¬ 
tainty, not only in the signal-depleted channels, but also in the most sensitive analysis channels for this 
search, (4 j, 3 b) and (4 j, 4 b). The channels with two f?-tags are used to constrain the leading uncertain¬ 
ties affecting the tt-flight-jets background prediction, while the channels with >5 jets and >3 f7-tags are 
sensitive to the uncertainties affecting the tt+HF background prediction. In particular, one of the main 
corrections applied by the fit is an increase of the tt + bb normalisation by approximately 20% relative to 
the nominal prediction by adjusting the corresponding nuisance parameter. This results in an improved 
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agreement between data and prediction in the (5 j, >4 b) and (>6 j, >4 b) channels, where the tt + bb 
process provides a significant, or, in the latter case, dominant background contribution.^ In addition, 
the corresponding uncertainty is reduced from the initial 50% down to about 16%. This correction is in 
agreement with that found in ref. [21]. However, in contrast with ref. [21], the tt + bb normalisation is not 
one of the leading uncertainties affecting this search, since the tt + bb background is subdominant in the 
(4 j, 3 b) channel, and it has a very different final discriminanf shape from fhaf of fhe signal in fhe (4 j, 
4 b) channel. 

As an illusfrafion, figure 7 provides a summary of fhe leading 15 sysfemafic uncerfainfies affecting fhe 
tt WbHc search, quantifying fheir impacf on fhe signal sfrengfh ju, bofh before and after fhe til, and 
displaying fhe constrainfs provided by fhe dafa on fhe associated nuisance paramefers. The pre-fil impacf 
on n is esfimafed by fixing fhe corresponding nuisance paramefer af Bq ± A0, where is fhe nominal value 
of fhe nuisance parameter and A0 is ifs pre-fil uncerlainly, and performing fhe fif again. The difference 
belween fhe defaull and modified //, A//, represenls fhe effecl of fhe sysfemafic uncertainly in quesfion on 
/r. The same procedure is followed lo eslimafe fhe posl-fit impacf on /r, bul fhe corresponding nuisance 
parameter is inslead fixed at 0 ± erg, where 8 is fhe filled value of fhe nuisance parameter and erg is ifs 
posl-fil uncerlainly. For reference, A/i = 0.05 corresponds lo ABR(f —> He) ^ 0.05%. 

Prior fo fhe til, fhe sysfemafic uncerfainfies wilh fhe largesl impact on are the leading uncertainty for 
light-jet tagging and the uncertainty on the tt background associated with the choice of parton shower 
and hadronisation models. The significant impact from light-jet tagging results from the large fraction of 
ff+light-jets background present in the (4 j, 4 b) channel, which peaks at high values of the final discrim- 
inanl, like fhe signal, and Ihus cannot be strongly constrained by the fit. Because of this, this uncertainty 
remains the leading one after the fit. In contrast, the uncertainty related to tt modelling is significantly 
constrained by the fit since it has a large impact (~5-16%) on the rt-i-light-jets background normalisation 
in the highly populated channels with two fj-tags. As a result, this uncertainty is ranked only fourth in 
importance after the fit, becoming comparable to uncertainties such as the choice of renormalisation scale 
for tt + bb, the leading uncertainty for c-jet tagging and the tt + cc normalisation. Of these, the nuisance 
parameter associated with the choice of the renormalisation scale for tt + bb is slightly pulled (by half 
of the prior uncertainty) to improve agreement with the data in the (4 j, 4 b) and (5 j, 4 b) channels. In 
these channels, this uncertainty causes variations of up to ~5% in the bin contents in some regions of the 
final discriminanf, i.e. disforfing its shape compared to that of the nominal prediction, but the sensitivity 
is not sufficient to constrain it significantly. The leading uncertainty from c-tagging causes small (few 
percent) distortions in the shape of the background, and also cannot be constrained by the fit. In contrast, 
the fit is sensitive to the tt + cc normalisation and the second-leading uncertainty for c-tagging,^ through 
the comparison of data and predictions across channels with different b-tag multiplicity, yielding results 
in agreement with the nominal predictions but with half the initial uncertainties. 

Other nuisance parameters have a smaller impact on the signal extraction and typically have small pulls 
or constraints. One exception is the nuisance parameter associated with the tt cross section, which affects 
the signal extraction indirectly through the existing small fraction of non-rt background, and after the fit 
is found to be consistent with the nominal prediction but is constrained owing to the large number of tt 
events. On the other hand, a slight pull is obtained for the nuisance parameter associated with one of 


® The overall change in tt + bb normalisation can be different across channels due to the different impact of other nuisance 
parameters affecting tt + bb modelling, such as that related to parton shower and hadronisation, which is changed by the fit by 
half of the prior uncertainty. 

^ The main effect of this uncertainty is a change in normalisation for the background with almost no effect on its shape in the 
final discriminant. 
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the uncertainties for the top quark pT and tt system px reweightings, which is used by the fit to improve 
agreement between data and prediction in the channels with two f?-tags but which has very small effect 
on the background prediction in the signal region. 

In the absence of a significant excess in data above the background expectation, 95% CL limits are set 
on BR(t ^ He) and BR(t ^ Hu). The observed (expected) 95% CL upper limits on the branching 
ratios are BR(t ^ He) < 0.56% (0.42%) and BR(t ^ Hu) < 0.61% (0.64%). These upper limits can 
be translated into corresponding observed (expected) limits on the couplings of \AtcH\ < 0.14 (0.12) and 
IT,<0.15(0.15). 


9.2. H -* yy 

A search for tt WbHq, H ^ yy published by the ATLAS Collaboration uses a data set corresponding 
to an integrated luminosity of 4.5 fb“^ at = 7 TeV and 20.3 fb“^ at = 8 TeV [17]. The event 
selection requires at least two reconstructed photon candidates and additional requirements to select t 
Wb decays. Events are categorised into two channels, leptonic and hadronic, depending on the W boson 
decay modes. The leptonic channel selects events with exactly one lepton {e or p), at least two jets, and 
at least one (j-tagged jet. The hadronic channel selects events with no reconstructed lepton, at least four 
jets, and at least one f^-tagged jet. In both channels, additional requirements are made to select events 
compatible with tt WbHq production by exploiting the invariant masses of the reconstructed top 
quark candidates. Finally, the diphoton mass (rriyy) distribution of the selected events is analysed using 
a sideband technique in order to estimate the background in the signal region, defined fo be 122 GeV < 
niyy <129 GeV. 

Based on fhe above sfrafegy, fhis search has essentially no discriminafion power befween tt —> WbHe 
and tt —> WbHu signals, because fheir selecfion accepfances are very close, alfhough nof identical. To 
facilifafe combining if wifh fhe ofher searches discussed in fhis paper, minor modifications fo fhe inpufs 
were made wifh respeef fo fhe published resulf, all having a negligible impaef on fhe resulf. They include 
updafes fo fhe tt cross-seefion uncerfainfy and fhe uncerfainfy model for Higgs branching rafios, as well as 
fhe separafe freafmenf of tt WbHe and tt —> WbHu signals faking info accounf fheir slighfly differenf 
accepfances. The besf-fif branching rafios obfained are BR(t —> He) = [0.22 + 0.26 (sfaf.) ± 0.10 (sysf.)]% 
and BR(t ^ Hu) - [0.23 ± 0.27 (sfaf.) + 0.10(syst.)]% under fhe assumptions fhaf BR(t Hu) - 0 
and BR(t ^ He) = 0 respectively. The observed (expeefed) 95% CL upper limifs on fhe branching 
rafios from ref. [17] remain fhe besf esfimafes, BR(f ^ Hq) < 0.79% (0.51%). The corresponding limifs 
on fhe couplings are |T,^//| < 0.17 (0.14). These limits can be understood as applying to the sum of the 
t He and t Hu decay modes, or only to one of them, if the other decay mode is assumed to have 
a branching ratio equal to zero. In the former case, BR(t ^ Hq) = BR(t ^ He) + BR(t ^ Hu) and 
\AtqH\ = + |T,„//P. 

9.3. H W+W-,T+T- 

The H WW* and H ^ tt decay modes are predicted to have significant branching ratios, of 21.5% 
and 6.3% respectively. The resulting signatures for signal events, tt WbHq —> W-W-W^bq and tt 
WbHq W-T-T^^bq, can be effectively exploited to suppress backgrounds in the case of multilepton 
final states. 
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(C) (d) 




(e) (f) 

Figure 4: tl —» WbHc, H ^ bh search; comparison between the data and prediction for the distribution of the D 
discriminant used in the (4 j, 2 b) channel (a) before the fit and (b) after the fit, in the (5 j, 2 b) channel (c) before 
the fit and (d) after the fit, and in the (>6 j, 2 b) channel (e) before the fit and (f) after the fit. The fit is performed 
on data under the signal-plus-background hypothesis. In the pre-fit distributions the tt —> WbHc signal (solid red) 
is normalised to BR(f ^ He) = 1% and the tt —> WbWb background is normalised to the SM prediction, while in 
the post-fit distributions both signal and tt —> WbWb background are normalised using the best-fit BR(f —> He). 
The small contributions from W/Z-i-jets, single top, diboson and multijet backgrounds are combined into a single 
background source referred to as “Non-fP’. The bottom panels display the ratios of data to either the SM background 
prediction before the fit (“Bkg”) or the total signal-plus-background prediction after the fit (“Pred”). The hashed 
area represents the total uncertainty on the background. 
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Figure 5: ff —» WbHc,H bb search; comparison between the data and prediction for the distribution of the D 
discriminant used in the (4 j, 3 b) channel (a) before the fit and (b) after the fit, in the (5 j, 3 b) channel (c) before 
the fit and (d) after the fit, and in the (>6 j, 3 b) channel (e) before the fit and (f) after the fit. The fit is performed 
on data under the signal-plus-background hypothesis. In the pre-fit distributions the tt —> WbHc signal (solid red) 
is normalised to BR(f ^ He) = 1% and the tt —> WbWb background is normalised to the SM prediction, while in 
the post-fit distributions both signal and tt —> WbWb background are normalised using the best-fit BR(f —> He). 
The small contributions from W/Z-i-jets, single top, diboson and multijet backgrounds are combined into a single 
background source referred to as “Non-fP’. The bottom panels display the ratios of data to either the SM background 
prediction before the fit (“Bkg”) or the total signal-plus-background prediction after the fit (“Pred”). The hashed 
area represents the total uncertainty on the background. 
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Figure 6: tl —» WbHc, H ^ bh search; comparison between the data and prediction for the distribution of the D 
discriminant used in the (4 j, 4 b) channel (a) before the fit and (b) after the fit, in the (5 j, >4 b) channel (c) before 
the fit and (d) after the fit, and in the (>6 j, >4 b) channel (e) before the fit and (f) after the fit. The fit is performed 
on data under the signal-plus-background hypothesis. In the pre-fit distributions the tt —> WbHc signal (solid red) 
is normalised to BR(f ^ He) = 1% and the tt —> WbWb background is normalised to the SM prediction, while in 
the post-fit distributions both signal and tt —> WbWb background are normalised using the best-fit BR(f —> He). 
The small contributions from W/Z+jets, single top, diboson and multijet backgrounds are combined into a single 
background source referred to as “Non-f?’. The bottom panels display the ratios of data to either the SM background 
prediction before the fit (“Bkg”) or the total signal-plus-background prediction after the fit (“Pred”). The hashed 
area represents the total uncertainty on the background. 
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Figure 1: tt ^ WbHc, H ^ bb search; the fitted values of the nuisance parameters for the most important sources 
of systematic uncertainty and their impact on the measured signal strength. The points, which are drawn conform¬ 
ing to the scale of the bottom axis, show the deviation of each of the fitted nuisance parameters, 6, from Oq, which is 
the nominal value of that nuisance parameter, in units of the pre-fit standard deviation Afl. The error bars show the 
post-fit uncertainties, erg, which are close to 1 if the data do not provide any further constraint on that uncertainty. 
Conversely, a value of erg much smaller than 1 indicates a significant reduction with respect to the original uncer¬ 
tainty. The nuisance parameters are sorted according to their post-fit effect on fj. (hashed blue area), conforming to 
the scale of the top axis, with those with the largest impact at the top. 
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Recently, the ATLAS Collaboration has published a search for ttH production with H WW*, tt and 
ZZ* [23] that exploits several multilepton signatures resulting from the leptonic decays of W and Z bosons 
and/or the presence of r leptons. That search considered five separate event categories depending on the 
number of reconstructed electrons or muons and hadronic t candidates, of which the following three are 
considered for a reinterpretation in the context of the tt WbHq search: 

• 2f’0Thad: two same-charge light leptons {e or ji) with no hadronic t candidates (rhad) and >4 jets 
with >1 /j-tagged jets. This channel is sensitive to the process tt WbHq —> i-i-qqqb2v. To 
further improve the sensitivity, this category is further subdivided into six subcategories depending 
on the flavour of fhe leptons and fhe number of jefs: {ee, jdfi, efj.)x{4 jefs, >5 jefs). 

• 3b: fhree lighf leptons wifh eifher >3 jefs of which >2 are /^-tagged, or >4 jefs of which >1 are 
fj-fagged. This channel is sensitive to fhe process tt —> WbHq —> i-b-i^qb3v. 

• 2/’lThad: two same-charge lighf leptons {e or ji), one Thad candidate, and >4 jefs wifh >1 /^-fagged 

jefs. This channel is sensitive to fhe process tt WbHq — > £-£-T'^qb3v. 

The fwo ofher cafegories considered in fhe ttH search in mulfilepfon final sfafes buf nof in fhis reinfer- 
prefafion are: 4£ (four lighf leptons wifh >2 jefs and >1 fj-lagged jefs) and 1^2Thad (one lighf lepton and 
fwo opposite-charge Thad candidates, wifh >3 jefs and >1 /7-lagged jefs), which have very small signal 
accepfance and/or poor signal-fo-background rafio due fo fhe large number of jefs misidenfified as Thad 
candidafes. The minimum number of jefs required in fhe 2f’0Thad cafegory is well mafched fo fhe number 
of parfons expecfed from tt WbHq W-W-W^bq signal evenfs, while fhe 3£ and 2/’lThad cafegor¬ 
ies effecfively require one or fwo jefs beyond leading order. Because of fhis, and fhe gain in branching 
rafio resulfing from only requiring fwo, as opposed fo fhree, W bosons decaying lepfonically, fhe 2(’0Thad 
cafegory dominates fhe sensifivify. 

The largesf background in fhe mosf sensifive cafegory, 2£4], is tt or single-fop producfion wifh one of 
fhe lepfons originating from a decay of a heavy-flavour hadron (“non-prompf lepfon”), followed by ttW 
producfion. Smaller confribufions arise from tt(Zly*), ttH, and diboson (primarily WZ) production, and 
dilepfon evenfs from tt and Zjy* wifh fhe wrong charge sign measured for one elecfron. In fhe remain¬ 
ing cafegories fhe non-prompf lepfon background decreases in imporfance, relafive fo fhe prompf-lepfon 
confribufions. The prompf confribufions are esfimafed using fhe simulafion, as fypically fhe relevanf 
processes (e.g., ttV, V - W,Z) have nof been measured af an accuracy exceeding fhaf of fheorefical pre¬ 
dictions. The non-prompf lepfon background predictions are compufed or validated using dafa confrol 
regions wifh similar lepton kinemafic selections buf wifh fewer jefs (fwo or fhree). Furfher defails are 
provided in ref. [23]. In fhe search for ttH, fhese regions have almosf no confaminafion from signal, buf 
fhe same is nof necessarily frue for fhe tt WbHq process, which is characferised by a lower jef mulfi- 
plicify fhan fhe ttH process. This could pofenfially lead fo an overesfimafe of fhe non-prompf background 
in fhis reinferprefafion, parficularly in fhe 2^0Thad cafegories, where fhe effecf is more pronounced. The 
3£ cafegory esfimafes fhe non-prompf lepfon background in fhe signal region by exfrapolafing, in ferms 
of lepfon qualify, from a dafa sideband region fhaf is sufficienfly deplefed in signal, using normalisation 
faclors obfained from fhe simulation. 

The non-prompt lepton background estimate for the 2f’0Thad categories is validated using simulation and 
control regions with two or three jets in which the lower px (subleading) lepton satisfies 10 GeV < px < 
30 GeV. Because fhe number of non-prompf lepfons increases significanfly as fhe lepfon px require- 
menf is lowered, fhis region is dominated by fhe non-prompf confribufion for any non-excluded value 
of BR(t ^ Hq). To exfracf fhe non-prompf background yields in fhese regions, binned likelihood fils 
fo fhe subleading lepfon px disfribulion are performed, accounting for non-prompf lepfon backgrounds. 
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Figure 8: ff —> WbHq, H WW*, tt search: comparison between the data and background prediction for the yields 
in each of the analysis channels considered before the fit to data. The expected tt WbHc and tt WbHu signals 
(dashed histograms) are shown separately normalised to BR(f —» Hq) - 1%. The sum of instrumental backgrounds 
originating from non-prompt leptons and lepton charge misidentification is denoted by “Non-prompt/q-misid”. The 
small contribution from rare processes such as tZ, ttWW, triboson, tttt and tH production are combined into a single 
background source denoted by “Rare”. The bottom panel displays the ratio of data to the SM background (“Bkg”) 
prediction. The hashed area represents the total uncertainty on the background. 


sources of prompt leptons (normalised to the SM expectation), and a potential tt WbHq contribution. 
The fitted tt —> WbHq contribution is in all cases found to be compatible with zero. The obtained best- 
fit non-prompt background yields in these low-pj control regions are used to normalise the simulation 
and obtain estimates for the background in the 2f’0Thad signal regions. These are found to be compatible 
with the nominal data-driven predictions within the stated uncertainties. The nominal predictions of the 
non-prompt backgrounds and their associated uncertainties are therefore used without modification in this 
reinterpretation. 

Figure 8 shows the event yields in each of the categories, which are used as input to the statistical analysis 
discussed in section 8. A table summarising the expected and observed yields can be found in appendix B. 
The best-fit branching ratio obtained is BR(t —> He) = [0.27 ± 0.18 (stat.) ± 0.21 (syst.)]% assuming that 
BR(f ^ Hu) = 0. A similar fit is performed for the tt WbHu search, yielding BR(t —> Hu) = [0.23 ± 
0.18 (stat.) ± 0.21 (syst.)]%, assuming that BR(f —> He) = 0. The observed (expected) 95% CL upper 
limits on the branching ratios are BR(t ^ He) < 0.79% (0.54%) and BR(t ^ Hu) < 0.78% (0.57%). The 
corresponding observed (expected) limits on the couplings are \AtcH\ < 0.17 (0.14) and \ AtuH\ < 0.17 (0.14) 
at 95% CL. 

9.4. Combination of searches 

The three searches discussed in the sections 9.1-9.3 are combined using the statistical analysis discussed 
in section 8. In this combination, the only systematic uncertainties taken to be fully correlated among 
the three searches are the tt cross section and the integrated luminosity for = 8 TeV data. The 
dominant uncertainties on the Higgs boson branching ratios primarily affect the tt WbHq, H ^ yy 
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and tt WbHq, H ^ bb searches. Other uncertainties such as those associated with leptons, jet energy 
scale and Zj-tagging should be partially correlated between the tt WbHq, H ^ yy search and the other 
two searches, but the differences in treatment between analyses (different lepton selection and px cuts, no 
uncertainty breakdown for jet energy scale and (^-tagging in the tt —> WbHq, H ^ yy search) make it 
difficult to account for correlations. However, given that the tt WbHq, H ^ yy search is completely 
dominated by the data statistics, the effect of this simplification is negligible. The uncertainties taken 
as correlated between the tt —> WbHq, H ^ bb and tt WbHq, H WW*,tt searches include 
those associated with lepton isolation, the leading (^-tagging and jet energy scale uncertainties, and the 
reweighting of top quark px and tt system px. The rest of the uncertainties are taken to be uncorrelated 
among the searches. This correlation scheme closely follows the procedure adopted in the combination 
of tlH searches by ATLAS [104]. 

The first set of combined results is obtained for each branching ratio separately, setting the other branching 
ratio to zero. The best-fit combined branching ratios are BR(t ^ He) = [0.22±0.10(stat.)±0.10(syst.)]% 
and BR(t ^ Hu) = [0.16 + 0.11 (stat.) ± 0.12 (syst.)]%. The difference between the central values of 
BR(f ^ He) and BR(f ^ Hu) originates from the ability of the H ^ bb search to probe both decay 
modes separately. A comparison of the best-fit branching ratios for the individual searches and their 
combination can be found in figure 9 for BR(f ^ He) and figure 10 for BR(t —> Hu). Figure 11 shows 
the CLs versus branching ratio for the combination. The observed (expected) 95% CL combined upper 
limits on the branching ratios are BR(t ^ He) < 0.46% (0.25%) and BR(t ^ Hu) < 0.45% (0.29%). 
The corresponding observed (expected) upper limits on the couplings are \AtcH\ <0.13 (0.10) and \AtuH\ < 
0.13 (0.10). A summary of the upper limits on the branching ratios obtained by the individual searches, 
as well as their combination, can be found in figure 12. 
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Figure 9: Summary of the best-fit BR(f —> He) for the individual searches as well as their combination, assuming 
that BR(f Hu) = 0. 


A similar set of results can be obtained by simultaneously probing both branching ratios. Although the 
tt WbHq, H ^ yy and tt WbHq, H WW*, tt searches are basically only sensitive to the sum of 
the two branching ratios, the tt WbHq, H ^ bb search has different sensitivity to each of them, and 
a simultaneous fit can be performed. The best-fit branching ratios obtained from the simultaneous fit are 
BR(f ^ He) = [0.34±0.22(stat.)±0.15 (syst.)]% andBR(t ^ Hu) = [-0.17±0.25 (stat.)±0.17 (syst.)]%. 
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Figure 10: Summary of the best-fit BR(f —> Hu) for the individual searches as well as their combination, assuming 
that BR(f —> He) - 0. 




Figure 11: (a) CLs versus BR(f He) and (b) CLs versus BR(f —» Hu) for the combination of the searches, 
assuming that the other branching ratio is zero. The observed CLs values (solid black lines) are compared to the 
expected (median) CLs values under the background-only hypothesis (dotted black lines). The surrounding shaded 
bands correspond to the 68% and 95% CL intervals around the expected CLs values, denoted by +lcr and +2cr, 
respectively. The solid red line at CLs =0.05 denotes the value below which the hypothesis is excluded at 95% CL. 


with a correlation coefficient of -0.84, as shown in figure 13. Figure 14(a) shows the 95% CL upper limits 
on the branching ratios in the BR(t ^ Hu) versus BR(t ^ He) plane. The corresponding upper limits on 
the couplings in the \AtuH\ versus plane can be found in figure 14(b). 
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Figure 12: 95% CL upper limits on (a) BR(f —» He) and (b) BR(f —> //m) for the individual searches as well as 
their combination, assuming that the other branching ratio is zero. The observed limits (solid lines) are compared to 
the expected (median) limits under the background-only hypothesis (dotted lines). The surrounding shaded bands 
correspond to the 68% and 95% CL intervals around the expected limits, denoted by ±lcr and +2cr, respectively. 
Because the asymptotic approximation is used in the calculation of CLs, the results of the H ^ yy search reported 
in this figure differ slightly from those published in ref. [17], which remain the most accurate results. 



Figure 13: Best-fit BR(f He) and BR(f —> Hu) and the corresponding 68% CL (solid) and 95% CL (dotted) 
regions for the combination of the searches. 


32 


















(a) (b) 

Figure 14: 95% CL upper limits (a) on the plane of BR(f ^ Hu) versus BR(f —> He) and (b) on the plane of |/1/„hI 
versus \A,ch\ for the combination of the searches. The observed limits (solid lines) are compared to the expected 
(median) limits under the background-only hypothesis (dotted lines). The surrounding shaded bands correspond to 
the 68% and 95% CL intervals around the expected limits, denoted by +lcr and +2cr, respectively. 


10. Conclusion 

A search for flavour-changing neutral current decays of a top quark to an up-type quark {q = u, c) and the 
Standard Model Higgs boson, where the Higgs boson decays to bb, is presented. The analysis searches 
for top quark pair events in which one top quark decays to Wb, with the W boson decaying leptonically, 
and the other top quark decays to Hq. The search is based on pp collisions at = 8 TeV recorded in 
2012 with the ATLAS detector at the CERN Large Hadron Collider and uses an integrated luminosity 
of 20.3 fb“^ Data are analysed in the lepton-plus-jets final state, characterised by an isolated electron 
or muon with moderately high transverse momentum and at least four jets. In this search the dominant 
background process is tt —> WbWb. To separate the signal from the background, the search exploits the 
high multiplicity of Z>-quark jets characteristic of signal events, and employs a likelihood discriminant 
that combines information from invariant mass distributions and the flavour of the jets. No significant 
excess of events above the background expectation is found, and observed (expected) 95% CL upper 
limits of 0.56% (0.42%) and 0.61% (0.64%) are derived for the t —> He and t —> Hu branching ratios 
respectively. 

Results from other ATLAS searches are also summarised, including a previous search probing H ^ yy 
decays, and a reinterpretation of a search for tJH production in multilepton final states, which exploits 
the H WW* ,TT decay modes. All searches discussed in this paper have comparable sensitivity, and 
thus their combination represents a significant improvement over the individual results. The observed 
(expected) 95% CL combined upper limits on the t He and t Hu branching ratios are 0.46% 
(0.25%) and 0.45% (0.29%) respectively. The corresponding observed (expected) upper limits on the 
\XtcH\ and \AtuH\ couplings are 0.13 (0.10) and 0.13 (0.10) respectively. Upper limits in the t He versus 
t Hu branching ratio plane, as well as best-fit branching ratios, are also reported. These are the most 
restrictive direct bounds on tqH {q - u, c) interactions measured so far. 
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Appendix 


A. Pre-fit and post-fit event yields in the tt WbHq, H bb search 

Table 3 presents the observed and predicted yields in each of the analysis channels for the tt —> WbHq, H —> 
bb search before the fit to data. Table 4 presents the observed and predicted yields in each of the analysis 
channels for the tt WbHc, H ^ bb search, after the fit to the data under the signal-plus-background 
hypothesis. 


B. Pre-fit event yields in the tt WbHq, H WW*, tt search 


Table 5 presents the observed and predicted yields in each of the analysis channels for the tt WbHq, H 
WW*, TT search before the fit to data. 
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4j,2b 

4j,3b 

4j,4b 

ti WbHc 

890 + 100 

394 + 54 

41.6 + 7.2 

tt WbHu 

851 +98 

339 + 49 

3.81 +0.71 


ff+light-jets 

77400 + 8100 

6170 + 860 

53 + 12 

tt + cc 

4900 + 2600 

680 + 370 

21 + 11 

tt + bb 

1870 + 990 

680 + 370 

44 + 23 

ttV 

121 +21 

15.5 + 2.9 

0.89 + 0.19 

tin 

30.5 + 4.2 

12.7 + 1.9 

1.91 +0.34 

fT+jets 

4700 + 1600 

217 + 78 

5.4 + 2.0 

Z+jets 

1080 + 450 

50 + 22 

0.90 + 0.50 

Single top 

4900 + 1400 

340 + 100 

6.8 + 2.3 

Diboson 

212 + 75 

11.5 + 4.1 

0.24 + 0.11 

Multijet 

1540 + 550 

100 + 36 

3.4+ 1.2 

Total background 

96800 + 9600 

8300+ 1100 

138 + 32 

Data 

98049 

8752 

161 



5j,2b 

5j,3b 

5j, >4b 

tt WbHc 

483 + 96 

242 + 50 

35.1 +7.7 

tt WbHu 

473 + 95 

217 + 46 

8.4 + 2.0 

ff+light-jets 

37600 + 6600 

3480 + 750 

61 + 18 

tt + cc 

4300 + 2300 

810 + 460 

43 + 28 

tt + bb 

1670 + 860 

890 + 470 

115 + 61 

tiV 

145 + 24 

26.5 + 4.5 

3.10 + 0.60 

ttH 

40.9 + 4.8 

22.3 + 2.9 

5.96 + 0.98 

VT+jets 

1850 + 790 

131+57 

5.8 + 2.7 

Z+jets 

400 + 200 

29+ 14 

1.47 + 0.76 

Single top 

1880 + 740 

195 + 78 

8.3 + 3.1 

Diboson 

96 + 41 

8.0 + 3.5 

0.40 + 0.19 

Multijet 

450 + 160 

68 + 24 

8.3 + 3.0 

Total background 

48400 + 7800 

5700+ 1100 

252 + 75 

Data 

49699 

6199 

286 



>6j,2b 

>6j,3b 

>6 j, >4 b 

tt WbHc 

267 + 68 

145 + 37 

31.1 + 8.3 

tt WbHu 

259 + 67 

132 + 34 

10.3 + 2.8 


ff+light-jets 

18800 + 4800 

2000 + 730 

52 + 40 

tt + cc 

3700 + 2000 

850 + 500 

79 + 46 

tt + bb 

1430 + 760 

970 + 520 

240+ 130 

ttV 

182 + 32 

44.6 + 8.1 

8.4+ 1.7 

ttH 

64.2 + 8.2 

39.8 + 5.4 

16.1 +2.6 

IT+jets 

880 + 440 

95 + 47 

8.5 + 4.5 

Z+jets 

180 + 100 

19+ 11 

1.5 + 0.9 

Single top 

840 + 410 

122 + 62 

11.9 + 6.2 

Diboson 

50 + 26 

6.0 + 3.0 

0.54 + 0.29 

Multijet 

176 + 62 

20.3 + 7.2 

0.93 + 0.50 

Total background 

26400 + 6100 

4200 + 1200 

420 + 160 

Data 

26185 

4701 

516 


Table 3: ff —> WbHq,H —> bh search; predicted and observed yields in each of the analysis channels considered. 
The prediction is shown before the fit to data. Also shown are the signal expectations for tt —> WbHc and tt —> 
WbHu assuming BR(f —> He) - 1% and BR(f —> Hu) - 1% respectively. The tt —> WbWb background is 
normalised to the SM prediction. The quoted uncertainties are the sum in quadrature of statistical and systematic 
uncertainties on the yields. 
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4j,2b 

4j,3b 

4j,4b 

tt WbHc 

155 + 11 

68.5 + 5.0 

7.23 + 0.56 


ff-H light-jets 

77300 + 1700 

6240 + 190 

56.1 +7.3 

tt + cc 

5600 + 1500 

810 + 210 

23.4 + 6.0 

tt + bb 

2420 + 360 

890 + 130 

54.1 + 7.5 

ttV 

122+ 19 

15.6 + 2.5 

0.90 + 0.15 

ttH 

30.9 + 3.6 

12.8+ 1.6 

1.92 + 0.26 

VT-i-jets 

4900+ 1100 

231 +55 

5.8 + 1.5 

Z-Hjets 

1040 + 390 

48 + 18 

0.78 + 0.32 

Single top 

5100+ 1000 

352 + 71 

7.1 + 1.5 

Diboson 

209 + 72 

11.6 + 4.0 

0.22 + 0.08 

Multijet 

1120 + 320 

75 + 22 

2.41 + 0.70 


Total 

98070 + 370 

8756 + 98 

159.9 + 7.4 

Data 

98049 

8752 

161 



5j,2b 

5j,3b 

5j,>4b 

tt WbHc 

85 + 10 

42.4 + 5.1 

6.17 + 0.77 


ff+light-jets 

37600 + 1200 

3570 + 160 

66.4 + 8.1 

tt + cc 

4700 + 1200 

970 + 240 

59+ 16 

tt + bb 

2140 + 310 

1150+ 160 

144 + 17 

ttV 

145 + 23 

26.6 + 4.2 

3.12 + 0.50 

ttH 

41.0 + 4.5 

22.4 + 2.6 

5.98 + 0.74 

fk+jets 

1940 + 560 

140 + 41 

6.0+ 1.8 

Z+jets 

380 + 170 

27.6 + 12.4 

1.38 + 0.63 

Single top 

2090 + 630 

219 + 66 

10.1 + 3.1 

Diboson 

93 + 39 

7.8 + 3.3 

0.37 + 0.16 

Multijet 

332 + 96 

46 + 13 

6.2+ 1.9 


Total 

49570 + 250 

6223 + 66 

308 + 11 

Data 

49699 

6199 

286 



>6j, 2 b 

>6j,3b 

>6 j, >4 b 

tt WbHc 

46.1 +4.2 

25.0 + 2.3 

5.41 + 0.52 


ff+light-jets 

18590 + 800 

2080 + 140 

54.2 + 8.4 

tt + cc 

3820 + 920 

980 + 240 

85+20 

tt + bb 

1860 + 270 

1260 + 170 

320 + 35 

ttV 

178 + 27 

43.7 + 6.8 

8.3 + 1.3 

ttH 

64.0 + 7.2 

39.6 + 4.5 

15.9 + 1.8 

IT+jets 

680 + 220 

75 + 24 

7.0 + 2.7 

Z+jets 

159 + 78 

16.8 + 8.3 

1.48 + 0.75 

Single top 

740 + 270 

108 + 40 

10.6 + 4.0 

Diboson 

48 + 23 

5.7 + 2.7 

0.51+0.25 

Multijet 

120 + 34 

13.9 + 4.0 

1.12 + 0.65 


Total 

26300 + 160 

4652 + 62 

508 + 22 

Data 

26185 

4701 

516 


Table 4: tt —> WbHc, H ^ bb search; predicted and observed yields in each of the analysis channels considered. 
The background prediction is shown after the fit to data under the signal-plus-background hypothesis. The quoted 
uncertainties are the sum in quadrature of statistical and systematic uncertainties on the yields, computed taking 
into account correlations among nuisance parameters and among processes. 
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eeA] 

ep4j 

A^A'4j 

tt WbHc 
ft WbHu 

^ 4 +0.S 
-0.6 

5.5 + 0.7 

15.0 

15.9 

11 1 +11 

97^ 

-i.n 


Non-prompt 

3.4+ 1.7 

12 + 4 

6.3 + 2.6 

q mis-id 

1.8+ 0.7 

1.4+ 0.6 

- 

ftW 

2.0 + 0.4 

6.2+ 1.0 

4.7 + 0.9 

tt(Zly*) 

0.75 + 0.20 

1.5 + 0.3 

0.80 + 0.22 

Diboson 

0.7 + 0.4 

1.9 + 1.0 

0.53 + 0.30 

tlH 

0.44 + 0.06 

1.16 + 0.14 

0.74 + 0.10 

Rare 

0 25 

-0.07 

0.72 + 0.05 

0 T,A +0 04 
'7-’^ -n ot 

Total background 

9.5+ 2.1 

25 + 5 

13.4 + 2.9 

Data 

9 

26 

20 



ee>5i 

e^i>5i 


tt WbHc 

2.9 + 0.6 

10.2 


tt WbHu 

2.7 + 0.8 

8 . 2 ;;® 

7 2 

-1.4 


Non-prompt 

2.3 + 1.2 

6.7 + 2.4 

2.9+ 1.4 

q mis-id 

1.1 +0.5 

0.85 + 0.35 


ttW 

1.4+ 0.4 

4.8 + 1.2 

3.8 + 0.9 

ftW) 

0.98 + 0.26 

2.1 + 0.5 

0.95 + 0.25 

Diboson 

0.47 + 0.29 

0.38 + 0.30 

0.7 + 0.4 

ttH 

0.73 + 0.14 

2.1 + 0.4 

1.41 +0.28 

Rare 

0.27 + 0.02 

0 79 +04)5 

^ - 0.04 

0.38 + 0.02 

Total background 

7.2+ 1.8 

17 + 3 

10.0 + 2.2 

Data 

10 

22 

11 


M 2f’lThad 

ft WbHc 

6.4 + 0.9 2.4 + 0.4 

tt WbHu 

5.2 

+ 1.0 2 1 

- 0.8 4..1 

+ 0.5 

- 0.4 


Non-prompt 

3.2 + 

0.7 

0 4 +')-6 

-0.4 

q mis-id 

- 


- 

tiW 

2.3 + 

0.7 

0.38 + 0.12 

ttiZIy*) 

3.9 + 

0.8 

0.37 + 0.08 

Diboson 

0.86 + 

0.55 

0.12 + 0.11 

ttH 

2.34 + 

0.35 

0.47 + 0.08 

Rare 

0 92 +0 0 ^ 
-0.06 

0 10 +“4)2 

Total background 13.7 + 2.3 

1.9+ 0.6 

Data 

18 

1 


Table 5: ti —> WbHq,H —> WW* ,tt search; predicted and observed yields in each of the event categories con¬ 
sidered. The prediction is shown before the fit to data. Also shown are the signal expectations for it —> WbHc and 
tt —> WbHu assuming BR(f —> He) =1% and BR(f —> Hu) - 1% respectively. The quoted uncertainties are the 
sum in quadrature of statistical and systematic uncertainties on the yields. 
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S. Istin^^®, J.M. Iturbe Ponce^"^, R. J. Ivarsson^^ W. Iwanski^^, H. Iwasaki^^, J.M. Izen'*^ 

V. Izzo^°^^, S. Jabbar^, B. Jackson*^^, M. Jackson^^, P. Jackson^ M.R. Jaekel^*^, V. Jain^, K. Jakobs^^, 

S. Jakobsen^*’, T. Jakoubek^^^, J. Jakubek^^^, D.O. Jamin^^^, D.K. Jana™, E. Jansen^^, R. Jansky^^, 

J. Janssen^^, M. Janus^^, G. Jarlskog*^ N. Javadov^^’^, T. Javurek^^, L. Jeanty^^, J. Jejelava^^®’“, 

G. -Y. Jeng^^*^, D. Jennens*^, P Jenni^^’", J. Jentzsch^^, C. Jeske^™, S. Jezequel^, H. J. Jia^^^, 

Y. Jiang^^*’, S. Jiggins’^, J. Jimenez Pena^^^, S. Jin^^®, A. Jinaru^^'’, O. Jinnouehi^^^, M.D. Joergensen^^, 
P Johansson^^^, K.A. Johns^, WJ. Johnson^^^, K. Jon-And^^^®’^^^^, G. Jones^™, R.W.L. Jones^^, 

T. J. Jones’^, J. Jongmanns^^^*, PM. JorgeK.D. Joshi^^, J. Jovieevie^^^^, X. 

A. Juste Rozas^^’^, M. Kaei^^’, A. Kaezmarska^^, M. Kado^^^, H. Kagan^^\ M. Kagan^^^, S.J. Kahn*^, 
E. Kajomovitz^^, C.W. Kalderon^™, S. Kama^°, A. Kamenshehikov^^*’, N. Kanaya^^^, S. Kaneti^^, 

V.A. Kantserov^^, J. Kanzaki^^, B. Kaplan^^°, E.S. Kaplan™^, A. Kapliy^^, D. Kar^‘^^‘^, K. Karakostas^*’, 
A. Karamaoun^, N. Karastathis™’^®^, M.J. Kareem^^, E. Karentzos^*^, M. Karnevskiy^^, S.N. Karpov^^, 

Z. M. Karpova^^, K. Karthik'^*’, V. Kartvelishvili^^, A.N. Karyukhin^^*^, K. Kasahara^^*’, E. Kashif^^^, 
R.D. Kass^^^ A. Kastanas^"^, Y Kataoka^^^, C. Kato^^^, A. Katre™, J. Katzy^^, K. Kawade^°^, 

K. Kawagoe’*’, T. Kawamoto^^^, G. Kawamura^^, S. Kazama^^^, V.E Kazanin^®^’^, R. Keeler^^^, 

R. Kehoe^°, J.S. Keller^^, J.J. Kempster’^, H. Keoshkerian^^, O. Kepka^^’, B.P Kersevan™, 

S. Kersten^^^, R.A. Keyes^^, E Khalil-zada^^, H. Khandanyan^^^^’^^^'’, A. Khanov^^^, 

A.G. Kharlamov^°^’^, T.J. Khoo^*, V. Khovanskiy^^, E. Khramov^^, J. Khubua^^*’’"', S. Kido^’, 

H. Y Kim^ S.H. Kim^^^, YK. Kim^i, N. Kimura^^^, o.M. Kind^^ B.T. King™, M. King^^^ 

S.B. King^^^, J. Kirk^^\ A.E. Kiryunin^*’^, T. Kishimoto^^, D. Kisielewska^^^, E. Kiss^^, K. Kiuehi^^°, 

O. Kivernyk^^^, E. Kladiva'™'’, M.H. Klein^^, M. Klein™, U. Klein™, K. Kleinkneeht^^, 

P. Klimek™®®’™^*’, A. Klimentov^^, R. Klingenberg^^, J.A. Klinger^^^, T. Klioutehnikova^**, 

E.-E. Kluge^^®, P. Kluit^®’, S. Kluth^‘’\ J. Knapik^^, E. Kneringer^^, E.B.E.G. Knoops^^, A. Knue^^, 

A. Kobayashi^^^, D. Kobayashi^^^, T. Kobayashi^^^, M. Kobel™, M. Koeian™^, P. Kodys^™, T. Koffas^^, 
E. Koffeman^*’^, E.A. Kogan^^°, S. Kohlmann^™, Z. Kohout^^^, T. Kohriki^^, T. Koi™^, H. Kolanoski^^, 
M. Kolb^*'’, I. Koletsou^, A.A. Komar^^’*, Y. Komori^^^, T. Kondo^^, N. Kondrashova^^, K. Kbneke^*, 
A.G. Kbnig^°^, T. Kono^^, R. Konoplieh'^®’^, N. Konstantinidis’^, R. Kopeliansky^^^, S. Koperny^^®, 

E. Kbpke^^, A.K. Kopp^^, K. KoreyP^, K. Kordas^^^, A. Korn’^, A.A. Korol^*’^’'^, I. Korolkov^^, 

E.V. Korolkova^^^, O. Kortner^*’^ S. Kortner^°\ T. Kosek^™, V.V. Kostyukhin^^, V.M. Kotov^^, 

A. Kotwal^^, A. Kourkoumeli-Charalampidi^^^, C. Kourkoumelis^, V. Kouskoura^^, A. Koutsman^^^^, 

R. Kowalewski^^^, T.Z. Kowalski^^®, W. Kozaneeki^^^, A.S. Kozhin^^°, V.A. Kramarenko™, 

G. Kramberger™, D. Krasnopevtsev^*, M.W. Krasny^°, A. Krasznahorkay^*’, J.K. Kraus^^, 

A. Kravehenko^^, S. Kreiss'^*^, M. Kretz^^'^, J. Kretzsehmar™, K. Kreutzfeldt^^, P. Krieger^^*, 

K. Krizka^^, K. Kroeninger^^, H. Kroha^*^^ J. Kroll^^^, J. Kroseberg^\ J. Krstie^^, U. Kruehonak^^, 

H. Kriiger^^, N. Krumnaek^^, A. Kruse^^^, M.C. Kruse^^, M. Kruskal^^, T. Kubota*^, H. Kueuk’^, 

S. Kuday4^ S. Kuehn^^ A. Kugel^^^ E Kuger^™, A. Kuhl^^^, T. Kuhl^^ V. Kukhtin^^ R. Kukla^^e, 

Y. Kulehitsky™, S. Kuleshov^^*’, M. Kuna*^^®’^^^*’, T. Kunigo^^, A. Kupeo^^^, H. Kurashige^’, 

YA. Kuroehkin^^, V. Kus^^’, E.S. Kuwertz^™, M. Kuze^^^, J. Kvita^^^, T. Kwan^™, 

D. Kyriazopoulos^^^, A. Ea Rosa^^^, J.E. Ea Rosa Navarro™‘^, E. Ea Rotonda^’^’^™, C. Eaeasta^^^, 

E. Eaeava^^^‘*’^^^’’, J. Eaeey™, H. Eaeker^^, D. Eaeour^**, V.R. Eaeuesta^^’, E. Eadygin^^, R. Eafaye^, 

B. Eaforge**^, T. Eagouri^^^, S. Eai^^, E. Eambourne^^, S. Eammers^^ C.E. Eampen’, W. EampF, 

E. Ean 9 on^^^, U. Eandgraf^*, M.P.J. Eandon’^, V.S. Eang^^®, J.C. Eange^^, A.J. Eankford^^^, E. Eanni^^, 
K. Eantzseh^\ A. Eanza^^'®, S. Eaplaee*°, C. Eapoire^*^, J.E. Eaporte^^^, T. Eari^^^*, 

E. Easagni Manghi™®’™’’, M. Eassnig^°, P. Eaurelli"^^, W. Eavrijsen^^, A.T. Eaw^^’, P. Eayeoek™, 

T. Eazovieh^’, O. Ee Dortz^*^, E. Ee Guirriee*^, E. Ee Menedeu^^, M. EeBlane^™, T. EeCompte®, 

E. Eedroit-Guillon^^, C.A. Eee™^^*, S.C. Eee^^^, E. Eee^, G. Eefebvre^*’, M. Eefebvre^™, E. Eegger^™, 

C. Eeggett^^, A. Eehan™, G. Eehmann Miotto^*^, X. Eei^, W.A. Eeight™, A. Eeisos^^^’^, A.G. Eeister^™, 
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M.A.L. Leite^^‘^, R. Leitner^^^, D. Lellouch^’^, B. Lemmer^^, KJ.C. Leney^^, T. Lenz^^, B. Lenzi^*^, 

R. Leone^, S. LeoneC. Leonidopoulos^^, S. Leontsinis^*^, C. Leroy^^, C.G. Lester^^, 

M. Levehenko^^^, J. Leveque^, D. Levin^^, LJ. Levinson^^^, M. Levy^^, A. Lewis^^°, A.M. Leyko^\ 

M. Leyton^i, B. h. h.L. Li^i, L. L. S. Li^^, X. Y. Z. Liang^^v^ 

H. Liao^"^, B. Liberti^^^^*, A. Liblong^^^, R Liehard^°, K. Lie^^^, J. Liebal^', W. Liebig^"^, C. Limbaeh^\ 

A. Limosani^^°, S.C. T.H. Lin*^, F. Linde^°^, B.E. Lindquist^^*, J.T. Linnemann^*^, 

E. Eipeles'^^, A. Eipniaeka^^, M. Eisovyi^**’, T.M. Eiss^^^, D. Eissauer^^, A. Eister^^^, A.M. Eitke^^^, 

B. Eiui5i.«^, D. Eiui^i, H. Eiu^'^, J. Eiu^^ J.B. Eiu^^^, K. Eiu*^ E. Eiu^^^ M. Eiu^^, M. Eiu^^b, y. Eiu^^b, 

M. Eivan^^^^’^^^'’, A. Eleres^^, J. Elorente Merino^^, S.E. Eloyd^^, E. Eo Sterzo'^\ E. Eobodzinska^^, 

P. Eoeh^, W.S. Eoekman^^^, E.K. Eoebinger^"*, A.E. Eoevsehall-Jensen^^, K.M. Eoew^^, A. Eoginov^’^, 

T. Eohse^^, K. Eohwasser^^, M. Eokajieek^^^, B.A. Eong^^, J.D. Eong^^^, R.E. Eong’^, K.A. Eooper^^\ 
E. Eopes^^^^, D. Eopez Mateos^^, B. Eopez Paredes^^^, I. Eopez Paz^^, J. Eorenz^®*’, 

N. Eorenzo Martinez^^ M. Eosada^^^, PJ. Ebsel^^*’, X. Eou^^®, A. Eounis^^^, J. Eove^, PA. Eove’^, 

H. Eu6°^ N. Eu^^ H.J. Eubattil3^ C. Eueii32ad32b^ ^ Eueotte^^ C. Euedtke^^ E. Euehring*’', 

W. Eukas^2^ L. Euminari'22a^ q Lundberg*'^^^’^^^'’, B. Eund-Jensen'^^^ P) Lynn^^, R. Eysak^^?^ 

E. Eytken^^, H. Ma^^, E.E. Ma22‘^, G. Maeearrone^^^ A. Maeehiolo^°\ C.M. Maedonald^^?^ g Macek’^, 
J. Maehado Miguens^ 22 , 1 26b ^ pj Maeina^*’, D. Madaffari*^, R. Madar^^, H.J. Maddoeks22, W.E Mader^^, 

A. Madsen'^^, J. Maeda^^^ S Maeland^^, T. Maeno^^, A. Maevskiy^^, E. Magradze^^, K. Mahboubi^^, 

J. Mahlstedt^*’^^ C Maiani^^e^ p; Maidantehik2^‘*, A.A. Maier^®^, T. Maier^°°, A. ]vlaio^ 26 a,i 26 b,i 26 d^ 

S. Majewski^^^, Y. Makida^^, N. Makovee^^’^ g Malaeseu^°, Pa. Maleeki^^, V.P Maleev^^^, E. Malek^^, 

U. Mallik^^, D. Malon^, C. Malone^^^, S. Maltezos^°, V.M. Malyshev^*’^, S. Malyukov^^, J. Mamuzie^2^ 
G. Maneini"^2^ g Mandelli^*’, E. Mandelli^^®, I. Mandic^^, R. Mandryseh^^^ j Maneira^^^^’^^^’’, 

E. Manhaes de Andrade Eilho^^'’, J. Manjarres Ramos^^^*’, A. Mann^®*’, A. Manousakis-Katsikakis^, 

B. Mansoulie^26^ g Mantifel*^^ M Mantoani^^, E. Mapelli^^, E. Mareh^^^‘^, G. Marehiori^*’, 

M. Mareisovsky^^2^ p p Marino^®^, M. Marjanovie^^^ P) g Marley^^, E. Marroquim^^®, S.P Marsden^^, 
Z. Marshall^^, E.E. Marti^’^ § Marti-Gareia^^’^ g Martin^*’, T.A. Martin^^o^ yj Martin^^, 

B. Martin dit Eatour^^, M. Martinez^^’^, S. Martin-Haugh^^\ V.S. Martoiu^^'’, A.C. Martyniuk’^, 

M. Marx^^^, E. Marzano^ 22 a^ ^ Marzin^*’, E. Masetti^^^ X Mashimo^^^, R. Mashinistov^^, J. Masik^^, 

A.E. Maslennikov^°^’^', I. Massa^°^’^°’’, E. Massa^°‘‘’^°'’, P. Mastrandrea^, A. Mastroberardino^ 2 a, 37 b^ 

T. Masubuehi^^^, P. Mattig^^s^ j Mattmann^^^ j Maurer2^'’, S.J. Maxfield^"^, D.A. Maximov^®^’^, 

R. Mazini^^\ S.M. Mazza^^^’^^'’, G. Me Goldriek'^^, S.P. Me Kee^^, A. MeCarn*^, R.E. MeCarthy^^^, 
T.G. MeCarthy2^, N.A. MeCubbin^^^, K.W. MeEarlane^^’*, J.A. Mefayden^*, G. Mehedlidze^^, 

S. J. MeMahon^^^, R.A. MePherson^^^’^, M. Medinnis^^, S. Meehan^^*, S. Mehlhase^*’®, A. Mehta^^, 

K. Meier^^^*, C. Meineek^°°, B. Meirose^\ B.R. Mellado Gareia'^^^^, E Meloni'^, A. Mengarelli^‘^‘*’^°'’, 

S. Menke^'^^, E. Meoni^^^, K.M. Mereurio^’^ § Mergelmeyer2i, P. Mermod^^, E. Merola'*’^®’^*’^*’, 

C. Meroni^'®, ES. Merritt^^, A. Messina^^^®’^^^*’, J. Metealfe^^, A.S. Mete^^^, C. Meyer^^, C. Meyer^^^, 
J-P. Meyer^26^ j Meyer^®^^ jj Meyer Zu Theenhausen^^®, R.P Middleton^^i^ S Miglioranzi^^^^’^^^‘^, 

E. Mijovic^\ G. Mikenberg^22^ ]y[ Mikestikova^^^^ M Mikuz’^, M. Milesi^*, A. Milie^°, D.W. Miller^^, 

C. Mills^^, A. Milov^22^ D.A. Milstead^^^‘*’'^^’’, A.A. Minaenko^^o^ y Minami^^^, I.A. Minashvili^^, 

A.I. Mineer'^**, B. Mindur^^^, M. Mineev^^, Y. Ming^23^ g ]y[ ]viir^2^ g p Mistry^^^, T. Mitani^’i^ 

J. Mitrevski^**®, V.A. Mitsou^^^^ A. Miueei^^, PS. Miyagawa^^^, J.U. Mjbrnmark^\ T. 

K. Moehizuki^^, S. Mohapatra^^, W. Mohr"^^, S. Molander^^^‘*’^^^’’, R. Moles-Valls^\ R. Monden^^, 

K. Mbnig^2^ p Monini^^, J. Monk^^, E. Monnier^^, A. Montalbano^^*, J. Montejo Berlingen^2^ 

E. Montieelli^^ S. Monzani'^^®’^^^*’, R.W. Moore^, N. Morange^'^^ P) Moreno^^^, M. Moreno Elaeer^^, 
P. Morettini^®^*, D. Mori^"^^, T. Mori^^^, M. Morii^^^ M Morinaga'^^, V. Morisbak^^^, S. Moritz*^, 

A.K. Morley^^°, G. Mornaeehi^*^, J.D. Morris’^, S.S. Mortensen^^, A. Morton^^^ g Morvaj^°2^ 

M. Mosidze^'*’, J. Moss^^^, K. Motohashi^^’^ g Mount^"^^, E. Mountrieha^^, S.V. Mouraviev^^’*, 
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EJ.W. Moyse^^, S. Muanza*^, R.D. Mudd'^, F. Mueller'°\ J. Mueller^^^, R.S.P. Mueller^°°, 

T. Mueller^^, D. Muenstermann^^, R Mullen^^, G.A. Mullier^^, FJ. Munoz Sanchez^^, 

J. A. Murillo Quijada^*, WJ. MurrayH. Musheghyan^"^, E. Musto^^^, A.G. Myagkov^^°’“‘^, 

M. Myska^^*, B.P. Nachman^^^, O. Nackenhorst^^, J. Nadal^^, K. Nagai^^*’, R. Nagai^^^, Y. Nagai*^, 

K. Nagano^^, A. Nagarkar^^\ Y. Nagasaka^^, K. Nagata^^°, M. Nagel^*’', E. Nagy^^, A.M. Nairz^°, 

Y Nakahama^*^, K. Nakamura^^, T. Nakamura^^^, I. Nakano*^^, H. Namasivayam^^, 

R. E Naranjo Garcia^^, R. Narayan^\ D.I. Narrias Villar^^^, T. Naumann^^, G. Navarro^^^, R. Nayyar^, 

H.A. Neal^^, PYu. Nechaeva^^, TJ. Neep^^, PD. A. Negri^^^‘‘’^^^’’, M. Negrini^*’®, 

S. Nektarijevic^°^, C. Nellist^^^, A. Nelson^^^, S. Nemecek^^^, P Nemethy^^*’, A.A. Nepomuceno^^^*, 

M. Nessi^*^’"^^, M.S. Neubauer^^^, M. Neumann^^^, R.M. Neves^^°, P Nevski^^, PR. Newman^^, 

D.H. Nguyen^, R.B. Nickerson^^*^, R. Nicolaidou^^^, B. Nicquevert^°, J. Nielsen^^’, N. Nikiforou^^, 

A. Nikiforov^^, V. Nikolaenko^^°’“^, I. Nikolic-Audit^°, K. Nikolopoulos^^, J.K. Nilsen'^^, P. Nilsson^^, 
Y. Ninomiya'^^, A. Nisati^^^®, R. Nisius^°\ T. Nobe'^^, M. Nomachi"^, I. Nomidis^^, T. Nooney^^, 

S. Norberg^^^, M. Nordberg^®, O. Novgorodova"^^, S. Nowak^°\ M. Nozaki^^, E. Nozka^^^, K. Ntekas^® 

G. Nunes Hanninger^*, T. Nunnemann^*’®, E. Nurse^^, F. Nuti^^, F. O’grady^, D.C. O’Neil^^^, 

V. O’Shea^^, EG. Oakham^^’^, H. Oberlack^®^, T. Obermann^\ J. Ocariz^*^, A. Ochi^^, I. Ochoa^^, 

J. P Ochoa-Ricoux 32 ^ S. Oda™, S. Odaka^^, H. Ogren^i, A. S.H. Oh^^, C.C. Obrni^, 

H. Ohmani*’^, H. Oide^^, W. Okamura^^^ H. Okawa^*^°, Y. Okumura^i, T. Okuyama^^, A. Olariu^^b, 
S.A. Olivares Pino^^, D. Oliveira Damazio^^, A. Olszewski^^, J. Olszowska^^, A. Onofre^^^^*’^^^®, 

K. Onogi^*’^, P.U.E. Onyisi^^’', C.J. Oram^^^®, M.J. Oreglia^\ Y. Oren'^^, D. Orestano^^^®’^^"^*’, 

N. Orlando^^^, C. Oropeza Barrera^^, R.S. Orr^^^, B. Oseulati^*^^’^*^*’, R. Ospanov^^, 

G. Otero y Garzon^^, H. Otono™, M. Ouehrif^^^‘*, F. Ould-Saada^^^, A. Ouraou^^^, K.P Oussoren**^’, 

Q. Ouyang^^^, A. Oveharova'^, M. Owen^^, R.E. Owen^^, V.E. Ozean^^^, N. Ozturk^, K. Paehal^^^, 

A. Paeheeo Pages^^, C. Padilla Aranda^^, M. Pagacova^^, S. Pagan Griso^^, E. Paganis^^^, F. Paige^^, 

P. Pais*^ K. Pajehelii^ G. Palaeinol 5 ‘’^ S. Palestini^o, M. Palka 3 ^^ D. Pallin^^, A. Palmai26ad26b^ 

YB. Pan^^^, E.St. Panagiotopoulou^°, C.E. Pandini**^, J.G. Panduro Vazquez^^, P. Pani^^^^’^^^'’, 

S. Panitkin^^, D. Pantea^^’’, E. Paolozzi^^, Th.D. Papadopoulou^°, K. Papageorgiou^^^, A. Paramonov^, 

D. Paredes Hernandez^^^, M.A. Parker^^, K.A. Parker^^^, F. Parodi^*’^*’^*’'’, J.A. Parsons^^, U. Parzefall^^, 

E. Pasqualueei^^^^, S. Passaggio^®^, F. Pastore^^^‘*’^^^'’’*, Fr. Pastore^’, G. Pasztor^^, S. Pataraia^^^, 

N.D. Patel^^°, J.R. Pater^^, T. Pauly^*’, J. Pearee^^^, B. Pearson^^^, E.E. Pedersen^^, M. Pedersen^^^, 

S. Pedraza Eopez^^’, R. PedroS.V. Peleganehuk'*’^’'^, D. Pelikan^^^, O. Pene^^^, C. Peng^^®, 

H. Peng^^'’, B. Penning^^, J. Penwell^^, D.V. Perepelitsa^^, E. Perez Codina^^^^^, 

M.T. Perez Gareia-Estan^^^, E. Perini^^‘‘’^^'’, H. Pernegger^*’, S. Perrella^*’^“’^°^’’, R. Pesehke^^, 

V.D. Peshekhonov^^, K. Peters^*’, R.F.Y. Peters^^, B.A. Petersen^*^, T.C. Petersen^^, E. Petit^^, 

A. Petridis^ C. Petridou^^^, P. Petroff^^^, E. Petrolo^^^^*, F. Petrueei^^^‘‘’^^^'’, N.E. Pettersson^^^, 

R. Pezoa^^*’, PW. Phillips^^', G. Piaequadio^^^, E. Pianori^™, A. Pieazio^^, E. Pieearo^^, 

M. Pieeinini^*^^'^*’’’, M.A. Piekering^^*’, R. Piegaia^^, D.T. Pignotti^^^, J.E. Pileher^^, A.D. Pilkington^^, 
A.W.J. Pin^ 4 ^ J. Pinai26a,i26b,i26d^ Pinamontii 6 ^“’^*^^‘=’"'^, J.E. Pinfold^, A. Pingel^^ S. Pires^°, 

H. Pirumov^^, M. Pitt^^^, C. Pizio^^®’^^'’, E. Plazak'^^®, M.-A. Pleier^^, V. Pleskot^^^, E. Plotnikova^^, 

P. Plueinskii 46 a,i 46 b^ ^ Pluth*’^, R. Poettgenl^ 6 ^’l'^'^^ E. Poggioli”^ D. Pohl^i, G. Poleselloi^i^ 

A. Poley^^, A. Polieieehio^’^’^^*’, R. Polifka^^^, A. Polini^°^, C.S. Pollard^^, V. Polyehronakos^^, 

K. Pommes^*^, E. Ponteeorvo^^^^, B.G. Pope^*’, G.A. Popeneeiu^^‘^, D.S. Popovie^^, A. Poppleton^°, 

S. Pospisil^^^, K. Potamianos^^, I.N. Potrap®^, C.J. Potter^^^, C.T. Potter^G. Poulard^*’, J. Poveda^*’, 
V. Pozdnyakov^^, M.E. Pozo Astigarraga^°, P. Pralavorio*^, A. Pranko^^, S. Prasad^*’, S. Prell^^, 

D. Priee^^, E.E. Priee^, M. Primavera^^^, S. Prinee^^, M. Proissl^^, K. Prokofiev^*’‘^, F. Prokoshin^^'’, 

E. Protopapadaki^^^, S. Protopopeseu^^, J. Proudfoot^, M. Przybyeien^^^, E. Ptaeek^^^, 

D. Puddu^^^‘*’^^^’’, E. Puesehel^^, D. Puldon^^^, M. Purohit^^’"^, P. Puzo^^^, J. Qian^^, G. Qin^^, Y. Qin^"^ 
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A. Quadt^"^, D.R. Quarrie^^, W.B. M. Queitsch-Maitland*^, D. Quilty^^, S. Raddum^^^, 

V. Radeka^^, V. Radescu^^, S.K. Radhakrishnan^^^, R Radloff^^^, P. Rados^^, F. Ragusa^^^’^^*’, 

G. Rahal'^^, S. Rajagopalan^^, M. Rammensee^*^, C. Rangel-Smith^^^, F. Rauscher^®*’, S. Rave^^, 

T. Ravenscroft^^, M. Raymond^*^, A.L. Read^^^, N.P. ReadiofF^^, D.M. Rebuzzi^^^®’^^^*’, 

A. Redelbach'^^, G. Redlinger^^, R. Reece^^’, K. Reeves^\ L. Rehnisch^^, J. Reichert^^^, H. Reisin^^, 

C. Rembser^*’, H. Ren^^^, A. Renaud^^^, M. Rescigno^^^^*, S. Resconi^^^, O.L. Rezanova^*’^’'^, 

P Reznicek^^^, R. Rezvani^^, R. Richter^®^, S. Richter^^, E. Richter-Was^^*’, O. Ricken^^, M. Ridel^*^, 

P Rieck^^, CJ. Riegel^^^, J. Rieger^"*, O. Rifki^^^, M. Rijssenbeek^^*, A. Rimoldi^^^®’^^^^, L. Rinaldi^*^^*, 

B. Ristic^^, E. Ritsch^*^, L Riu^^, E. Rizatdinova^^^, E. Rizvi’^, S.H. Robertson^^’^, 

A. Robichaud-Veronneau^^, D. Robinson^^, J.E.M. Robinson^^, A. Robson^^, C. Roda^^^^’^^^'’, S. Roe^*^, 

O. R 0 hne^^^, A. Romaniouk^^, M. Romano^^^’^'^*’, S.M. Romano Saez^^, E. Romero Adam^^’, 

N. Rompotis^^^, M. Ronzani^*, E. Roos^°, E. Ros^^’, S. Rosati^^^^*, K. Rosbach^^, P Rose^^’, 

PE. Rosendahl^^, O. Rosenthal'^^, V. Rossetti^^^‘‘’^^^'’, E. Rossi^*^^^’^°^'’, E.P Rossi^*^®, J.H.N. Rosten^*, 
R. Rosten'^^, M. Rotaru^^*’, E Roth^^^, J. Rothberg^^^, D. Rousseau'C.R. Royon'^^, A. Rozanov^^, 

Y. Rozen'^^, X. Ruan'^^‘^, E. Rubbo'^^, E Rubinskiy^^, V.E Rud^^, C. Rudolph^^, M.S. Rudolph'^*, 

E. Ruhr"'^, A. Ruiz-Martinez^°, Z. Rurikova"'^, N.A. Rusakovieh*’^, A. Rusehke'°°, H.E. Russell'^^, 
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